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Abstract 

In this article we review a less known unperturbative and powerful many-body 
method in the framework of classical statistical mechanics and then we show how it 
works by means of explicit calculations for a nontrivial classical model. The formalism 
of two-time Green functions in classical statistical mechanics is presented in a form 
parallel to the well known quantum counterpart, focusing on the spectral properties 
which involve the important concept of spectral density. Furthermore, the general in- 
gredients of the classical spectral density method (CSDM) are presented with insights 
for systematic nonperturbative approximations to study conveniently the macroscopic 
properties of a wide variety of classical many-body systems also involving phase tran- 
sitions. The method is implemented by means of key ideas for exploring the spectrum 
of elementary excitations and the damping effects within a unified formalism. Then, 
the effectiveness of the CSDM is tested with explicit calculations for the classical 
(/-dimensional spin-S' Heisenberg ferromagnetic model with long-range exchange in- 
teractions decaying as r^P (p > d) with distance r between spins and in the presence 
of an external magnetic field. The analysis of the thermodynamic and critical proper- 
ties, performed by means of the CSDM to the lowest order of approximation, shows 
clearly that nontrivial results can be obtained in a relatively simple manner already to 
this lower stage. The basic spectral density equations for the next higher order level 
are also presented and the damping of elementary spin excitations in the low tempera- 
ture regime is studied. The results appear in reasonable agreement with available exact 
ones and Monte Carlo simulations and this supports the CSDM as a promising method 
of investigation in classical many-body theory. 
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1. Introduction 

At the present a wide variety of methods exists to calculate the macroscopic equilibrium 
and nonequilibrium quantities of many-body systems. However, their potentialities and ef- 
ficiencies differ sensibly especially when small parameters are absent for intrinsic reasons 
and the ordinary perturbative expansions appear inadequate. There exist, for instance, many 
systems of experimental and theoretical interest which involve strongly interacting micro- 
scopic degrees of freedom and exhibit anomalous behavior of the specific heat, suscepti- 
bility and other macroscopic quantities which, in turn, indicate various phase transitions 
to occur in these systems. Therefore, there is always a great interest to search for reliable 
methods going beyond the conventional perturbation theory to describe correctly nonper- 
turbative phenomena and, especially, the various anomalies in the thermodynamic behavior 
of systems near a phase transition. In this connection, the two-time Green function (GF) 
technique constitutes one of the most powerful tools in quantum statistical mechanics and in 
condensed matter physics to explore the thermodynamic and transport properties of a wide 
variety of many-body systems. Within this framework, the equations of motion method 
(EMM) and the spectral density method (SDM) allow to obtain reliable approximations to 
treat typically nonperturbative problems LL^^^^. 

When one mentions the two time GF technique, one refers usually to their quantum 
many-body version which is well known from long time and widely and successfully ap- 
plied in quantum statistical physics. Nevertheless, the pioneering introduction of the two- 
time GF's and the EMM in classical statistical mechanics by Bogoljubov and Sadovnikov 
||6 1 has opened the concrete possibility to describe classical and quantum many -body sys- 
tems on the same footing. In this context, it is worth mentioning that, in many physical 
situations (when the quantum effects are negligible), the use of the classical two-time GF 
formalism may offer substantial advantages especially from the computational point of view 
because in the calculations one handles only functions and not operators. In any case, the 
strong impact of the two-time GF technique in many-body physics and the continuous ef- 
forts for obtaining better approximations are well known and, hence, it is unnecessary to 
stress again its effectiveness and potentialities. Rather, the less known SDM merits further 
remarks. 

The SDM was formulated some decades ago by Kalashnikov and Fradkin for quantum 
many-body systems |4||5|. It is sufficiently simple and has been apphed to a wide variety of 
quantum systems, also exhibiting phase transitions |5|, as superconductors fllEl, magnetic 
BjiSj^TJ and bosonic |8 1 systems, strongly correlated electron compounds t5i|9in0l^D> and 
so on. It appears more effective than the EMM for a direct calculation of two-time GF's, 
not only because it allows to obtain reliable and systematic nonperturbative approximations, 
but also because it assures the validity of a number of sum rules irrespective either to the 
interaction strength or any other characteristics inherent to the system under study. 

Next, a classical version of the SDM (CSDM) has been also formulated in extensive 
|[T2l fT3l IT4I and nonextensive fTSll classical statistical mechanics and applied to classical 
magnetic systems with short and long-range exchange interactions. The last method, strictly 
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related to the classical two-time GF technique fE\, offers a robust technical instrument for 
systematic and well tested nonperturbative approximations to explore the equilibrium and 
transport properties of classical many-body systems as well, parallel to the quantum coun- 
terpart. Unfortunately, both the mentioned methods in the classical context have not re- 
ceived the due consideration and further developments and applications are desirable. 

In this article, also with a pedagogical aim in mind, we first review the basic formalism 
of two-time GF's in classical statistical mechanics in a form strictly parallel to the quantum 
better known counterpart emphasizing the concept of spectral density (SD), quite important 
for next purposes, and its relation to the two-time GF's. Besides, as a convenient alternative 
to the EMM for calculation of the GF's, the basic ingredients of the CSDM aie presented in 
some details, together with useful insights for systematic nonperturbative approximations. 
Then, we clarify how this method can be conveniently used to explore macroscopic prop- 
erties of classical many-body systems. We introduce also the key ideas for studying, in a 
conceptually transparent way, the elementary excitation spectrum, or dispersion relation, 
and the damping effects in a unified way within the formalism of CSDM. The general pre- 
sentation of the method in the classical many -body framework will be performed in a form 
strictly related to the well known quantum counterpart U|2l|5l, giving the possibility to use 
powerful calculation techniques, approximation methods and the basic terminology which 
are typical of the quantum many-body theory. 

Next, we will show, with explicit calculations, how the CSDM works emphasizing its 
effectiveness already to the lowest order of approximation. As a nontrivial theoretical lab- 
oratory, we choose the classical t/-dimensional spin-S Heisenberg ferromagnet with long- 
range interactions (LRI's), decaying as r^^ {p > d) with distance r between spins, in the 
presence of an external magnetic field. For this model, besides some partial exact results 
llT6llT7lfT8llT9ll20l . a lot of new information about thermodynamic and critical properties 
have been obtained by means of Monte Carlo (MC) simulations |2r "2T|. As we shall see, 
the CSDM allows us to obtain, in a relatively simple way, untrivial results which appear to 
be in a reasonable agreement with the available exact ones and MC simulations supporting 
the SDM as a promising method of investigation also in classical statistical mechanics. 

The paper is organized as follows. In Sec. 2, we introduce the two-time GF's and the 
related SD's in classical statistical mechanics for two arbitrary dynamic variables. Section 3 
is devoted to the methods of calculation focusing on the general formulation of the CSDM. 
Here, the method is properly implemented for exploring also the damping of the classical 
excitations within the same basic formalism. The method at work is presented in Sec. 4 
with a study of the thermodynamic properties and the critical phenomena of the classical 
cf-dimensional spin-5 Heisenberg ferromagnetic model with LRI's. Here, after a summary 
of known theoretical and simulational information about the model for a next checking of 
the effectiveness of the CSDM, we present explicit calculations performed by using the 
general many-body formalism developed before. As a successive useful comparison, we 
first formulate the problem in the context of the EMM for the two-time GF's appropriate 
to the spin model under study and the basic equations in the Tyablikov-like approximation 
(classical version of the well known quantum Tyablikov one |2|) are presented. Then, the 
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basic SD and related moment equations (ME's) are derived, the thermodynamic and critical 
properties of the model are presented working in the so called one-pole approximation for 
the SD and a comparison with the Tyablikov-like results is performed. To underline the 
effectiveness of this lowest order approximation and the nature of difficulties arising to 
higher orders, we deduce the basic equations for the two-poles approximation. As a further 
insight, we use the implemented CSDM to obtain, in a unified and relatively simple way, 
the dispersion relation and the damping of classical spin excitations. Explicit results are 
given (for the short-range FM chain) in the low-temperature limit. Finally, in Sec. 5 some 
conclusions are drawn. 



2. Two-Time Green Functions and Spectral Densities 
in Classical Statistical Mechanics 

2.1 Basic Definitions and Spectral Representations 

In strict analogy with the quantum case (2|2l|3l|23l, we define |6l| the two-time retarded 
(v = r) and advanced (v = a) GF's in classical statistical mechanics for two arbitrary dy- 
namical variables A and B as' 

G^^^{t,t') =Q,{t -t'){{A{t),B{t')}) ^ {{A{ty,B{t')))^ , {v = r,a), (2.1) 

where Qr{t - 1') = 6(f - 1'), Qa{t - 1') = -6(f' - 1) and 6(x) is the usual step function. 
In Eq. (I2.lt . (...) denotes an equilibrium or a nonequilibrium ensemble average |23| and 
{A,B} is the Poisson bracket of A and B. In the following, we refer only to the equilibrium 
(generalized) ensemble average. In this case one can easily prove 1231 that the two-time 
GF's (I2.lt and the corresponding correlations functions depend on times t, t' only through 
the difference t — t' . So, one can write 

gJ2(?-0 = {{A{t-t');B))^ = {{A;B{t'-t))),. (2.2) 

In the previous equations, the dynamical variables A and B depend on time via the conjugate 
canonical coordinates {q{t),p{t)) = {qi{t),..,q^{t);pi{t),..,p^{t)), (5\£ is the number of 
degrees of freedom of the classical system under study), X{t) = e'^X{0) withX =A,B,q,p, 
L = i{H,...} is the Liouville operator^ and H is the Hamiltonian. Here, e'^' acts as a classi- 
cal time-evolution operator which transforms the dynamical variable X{0) =X{q{0),p{0)) 
at the initial time t = into the dynamical variable X(f) =X{q{t),p{t)) at the arbitrary time 
t. Of course, the time evolution of the generic dynamical variable X{t) is governed by the 
well known Liouville equation of motion (EM) 

_^ = {X(0,//}. (2.3) 

'The classical retarded GF, as in the quantum one, has a direct physical meaning which arises from the 
theory of linear response for transport phenomena in classical statistical mechanics. The advanced one is 
introduced only for a full development of the two-time GF formalism. 

^In literature [23] one finds also the alternative definition L = —i{H,...}. However, the final conclusions 
remain unchanged. 
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With the previous definitions, one can show |6l| that the two-time correlation function 
FAB{t,t') = pAsit - t') = {A{t)B{t')) = (A(x)B) = (AB(-t)), with x = f - is related to 
the classical GF's (I2.2t by the relation 

dx 



G^^>{x) = Pev(x)3; (A(T)B) = pev(T) {{Aiz),H}B) , (2.4) 



where p = (KbT) \ T is the temperature and Kb is the Boltzmann constant. In particular, 
we have also 

({A(x),B}) = (A(T)B) = P {{A{z),H}B) , (2.5) 

which connects the Poisson bracket of two dynamical variables and the corresponding dy- 
namical correlation function. 

Working, as assumed before, within equilibrium Gibbs ensembles, one can introduce 
for (x) and Fab(t) the Fourier transforms 



FABiT) = ^PAB{(0)e-'''\ (2.7) 

where G^^ (co) = ( (A (x) ; B) )y ^ will be named the v-GF's of A and B in the co-representation 
and 7^AB(tt>) = (A(x)B)^ is called the classical spectral intensity of the time-dependent corre- 
lation function Fab(x), with 3((o) = J^^ dTe'""3{x). Then, using Eq. (I2.4t and the integral 
representations 

r+°° dx e~'^ r+°° dx ■ 

Q(^) = / ' ^ ^ 0+; 5W = / ^e'^, (2.8) 

for the step function and the Dirac 5-function, the Fourier transforms of the two-time GF 
(12.11) for two dynamical variables A and B can be expressed in terms of the corresponding 
spectral intensity as: 



GabH-I ^a)_co'+(_i)v,-e ' ' (2-9) 



where the symbol (—1)^ means +1 if v = r and — 1 if v = a. It is interesting to compare 
Eq. (12.91) with the quantum corresponding one for two operators A and B fni2ll3l l23l 



(V).. f^-dco' {l+^e-f'^^')FAB{oy') ^ „+ 



where r] = — 1 and r] = +1 for definition of quantum two-time GF's with commutator or 
anticommutator, respectively, and H is the reduced Planck constant. Notice that formally. 
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as expected for internal consistency, the function C((o) = Pco or Q^{oi) = (1 + r\e^^^)/h 
characterizes the classical or quantum nature of the problem under study, respectively. 

In analogy to the quantum case |3|5l, we now introduce the time-dependent classical 
spectral density (CSD) for A and B fT2lfT3lfT4ll24l 



A^fi(x) = /({A(x),B}). 



Taking its Fourier transform, we get 



(2.11) 



Aab((o) = /({A(x),B})„ = / dTe'^^'AABiT) = P(0^ab((o). 

J — oo 

Hence, from Eq. (I2.9t . one immediately obtains the spectral representation 

'■+~ doi' Aab((o') 



(2.12) 



J — < 



271 (0-(0'+(-l)^/£ 



0^ 



(2.13) 



for the two-time GF's ( 12.11 ) in terms of the corresponding CSD Aab(c)) in the (O-representation. 
Also the dynamical correlation function (A(t)B) can be easily expressed in terms of Aab{(o). 
From Eqs. i2.5l . ( 12.111 ) and ( 12.121 ). we obtain indeed 



{Aim 



271 Pro 



(2.14) 



Eqs. (I2.13t and (I2.14t assume a particular importance for our purposes. If one determines 
independently the Aab((o), these equations allow us to obtain the time-dependent correla- 
tion and Green functions and therefore the macroscopic properties of classical many-body 
systems. We will show below that, also in classical statistical mechanics, it is possible to 
construct a formalism which allows for a systematic calculation of the CSD. 

From Eqs. (12.1 U - dTTH i some immediate exact results can be easily obtained. First, 
Eqs. (I2.11t and (I2.12t yield 



Besides, from (12. 14t . it follows 



'<iro, 

^Aab((o) 



' <iro Aab(co) 
271 Pro 



:/({A,B}) 



(2.15) 



(AB) 



(2.16) 



The relations (I2.15t and (I2.16t are formally exact and constitute useful examples of the so 
called sum rules of the CSD, AAB(ro), which have great relevance for physical consistency 
of practical calculations and approximations. Combining now Eqs. (12. 13t and (I2.15t . one 
can easily prove another general result which plays an important role for calculation of the 
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GF's. As CO ^ oo we have indeed 

'■+~<ico' Aab((o') 



CO 



-1 



271 ^ a)'-(-l)v/E 



1 f+°°d(0 , , / 
- / — Aab((0 

C0 7-oo 271 



co'-(-l)^/8 1 ((o'-i-lYie 

1 + ^ — — + ^- — ^ + - 

CO 2 



C02 



i{{A,B}) 1 /■+- Jco' . ,w / . ^v N /I 
^' + — / — Aab(co') (co'-(-1)^/8) + o(-^ 



CO 



co^ J-oo 2% 



So, we can write 



gS(c«) 



CO 



at 

-a 



CO 



-1 



(a>2) 



if 
if 



{{A,B})=0, 



(2.17) 



(2.18) 



which provide a relevant boundary conditions for the v-GF's. 

Let's come back now to the relations (12.131) for classical retarded and advanced GF's 
in the co-representation. As in the quantum counterpart [1, 2|, using the definitions (12.11) 
and (12. 2t . one can easily prove that G^g(co) and G^'^(co), analytically continued in the 
co-complex plane, are analytical functions in the upper and lower half-plane, respectively. 
Then, combining these two analytical functions, one can construct a single function 
Gab(g>) = Jl^^dte'"^'GABit) of complex CO such that 



Gab((o) 



Im CO > 
Im CO < 0. 



(2.19) 



Hence, Eq. (I2.13t yields for Gab (co) the spectral representation 

Jco'Aab(co') 



Gab (co) 



271 CO - co' 



(2.20) 



This function is analytical in the whole complex co-plane with a cut along the real axis 
where singularities for Gab{c^) may occur. It is worth noting that, in terms of Aab(g>), 
no formal differences exist for the spectral representations of G^^(co) and GabCo) in the 
classical and quantum context. Hence, all the developments already known in quantum 
many-body theory remain formally valid for the classical case. So, one has immediately the 
important exact relation 



AAB(ro) = / [Gab((0 + i£) - Gab((0 - : 



(2.21) 



which express the CSD in terms of the related two-time GF's in the co-representation. This 
allows us to state also that the cut for Ga6((o) along the real axis in co-complex plane is 
determined by Eq. (12.211) and its singularities are the points of the real axis where the 
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condition A^g((o) 7^ is satisfied. For the spectral intensity of classical systems, Eqs. (I2.12t 

. Gab (co + ie) - Gab (co - ie) 



and (12.2 It yield 



Fab((o) = (A(x)B)„ = /- 



p(0 



(2.22) 



Of course, other known quantum relations are formally valid for classical many-body the- 
ory. For instance, when AAB(tt)) is real, the classical &amer-Kronig relations (classical 
dispersion relations) between the real and imaginary parts of G^g((o) and G^g(a)) for real 
values of co are true 



4)(v) 



71 



,Im_G^ 
CO' -co 



d(o' 



(2.23) 



where the symbol p denotes the main part of the integral defined as 



X — Xq 8^0+ 



X — Xq 



xo+^ X — Xq 



We have also 



AAB(co) = -2(-irimGj2(co) 



and, in particular. 



Aab(co) 



-2ImGS((0). 



(2.24) 



(2.25) 



(2.26) 



Many physical information about the macroscopic properties of classical many body 
systems can be obtained from the analytical properties of the GF's defined above. In the 
Subsec. 2.2 we will derive a spectral decomposition for Aab(co), and hence for G^^^(co) and 
F^b(co), from which it becomes possible to analyze the nature of the classical GF singular- 
ities and their connection with the physical quantities. 



2.2 Spectral Decompositions, Classical Oscillations and Damping 

In this section we will show that, in strict analogy with quantum case, also in classical 
statistical mechanics an exact spectral decomposition for the CSD is true. With this aim, 
following Refs. [25 26 1, we introduce the Hilbert space of the classical dynamical variables 
with a scalar product defined as 

{A\B)'^= Z{A*B) = J dre-^"A*{q,p)B{q,p), (2.27) 

where dV = dqdp and Z = J dFe^^^ is the (generalized) partition function for the prob- 
lem under study. In this space one can consider the eigenvalue equation for the Hermitian 
operator L 

L|\^^) =co^|\^^), (2.28) 
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or, equivalently, in terms of the wave functions \\fic{q,p) corresponding to the vectors 

L\^k{q,p) =(iik^k{q,p)- (2.29) 



In Eqs. (I2.28t and (I2.29I) . \\\fk) (}\ik) and (Ot denote the eigenvectors (eigenfunctions) and 
the eigenvalues of the Liouville operator L. 

An important property of the Liouville operator is that all its eigenfunctions are (in 
general) complex 1 25 1 and one can promptly prove that, if is an eigenfunction of L with 
eigenvalue cO/t, then \|/| is also an eigenfunction of L with eigenvalue —00^;. 

At this stage, if we assume that {\|/yt} is a complete set of orthonormal eigenfunctions 
i{\\fkWk') = ^ {^li^j p)^k{Q, p)) = 5i:,i:'), for two arbitrary dynamical Variables A (g,p) and 
B{q,p) we can consider the expansions: 



KP^q) = Zk{^*k\A)^*k{q,p), 
BiP,q)='LkWk\B)m{q,p)- 

On the other hand, from the relations (I2.5t . (I2.11t and (I2.12t . it follows that 

/+00 
jTe'""(A(x)B) 
-oo 



(2.30) 



(2.31) 



Then, taking into account the definition (I2.27t . the expansions (I2.30t and the orthonormality 

(2.32) 



and completeness conditions in the Hilbert space, Eq. (I2.31t yields 

Aab («) = 27ip(0Z- 1 £ (\^| |A ) (\^|. |B ) 6((0 - ©0 , 



where we have used the relations 



and 

^^^Km-m,)x^ 27i5((0 -(Ok). (2.34) 

-oo 

Equation (I2.32t is the desired spectral decomposition, or classical Lehmann representation, 
for the Aab((o) to be compared with the known quantum analog [4 ^. Of course, if B = A*, 
we have 

Aaa' («) = 27ip(0Z- ' ^ I (^^^ |A ) |2 6 ((0 - (0,) , (2.35) 

k 

which is a real quantity. Equations (I2.32t and (I2.35t express formally Aab((o) as an infinite 
sum of appropriately weighted 6-functions. It may happen that, for special physical cases, 
only a finite number of terms in the sum of Eqs. (I2.32t and (I2.35t are different from zero. 
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In such cases, the many-body problem may be solvable. Usually, however, Aab((o) will be 
a continuous function of real frequency co. 

It is easy to see that Aab((o) can be also written in the compact form 

Aab(co) =27ip(o(B6((0 + L)A) =27i/({5((0 + L)A,B}), (2.36) 



which, in view of definition (12.271) . is just the operatorial version of the spectral decompo- 
sition (I2.32t. 



An important consequence of the general expression (12.321) for Aab((o) is that the 



GAB(ro)> defined by Eqs. ( 12. 19t and (I2.20t . can be written as 

Gab((0) = ((A(x);B))„ = PZ-1£(V, \B) {Y^ |A) (2.37) 

Hence, the eigenvalues of the Liouville operator L = i{H, ...}, which is related to the energy 
of the system, are just the real poles of the Green function Gab((o). On the other hand, from 
Eq. (I2.14t . one has for the x-dependent correlation function Fab{'^) 

Fab(x) = (A(T)B) =Z-'£(v,|B)(v^|A)e-'»*^ (2.38) 



Then, in strict analogy to the quantum case, Eqs. (I2.37t and (I2.38t suggest that the real 
poles of Gab(g>)> the eigenvalues (Ok of the Liouville operator, represent the classical 
frequency spectrum of undamped oscillations (elementary or collective depending on the 
physical nature of the dynamical variables A and B) in the system. A characteristic fre- 
quency of this type will be called the oscillation dispersion relation. Thus, in case of polar 
singularities, the real poles of Gab(g>) represent undamped oscillations of the system asso- 
ciated to the time-dependent correlation function for the dynamical variables A and B. It 
is worth noting that, in general, the 8-poles in the exact classical Lehmann representation 
(12.3 It of the spectral density are expected lying infinitesimally close, therewith defining a 
continuous function AAB(tt>)- It is then obvious that the damped oscillation concept works 
only under the basic presumption that the CSD exhibits some pronounced peaks, which 
have to be identified with oscillation poles for Gab{c^)- The widths of these peaks are a 
direct measure of damping and, therefore, of the finite lifetime of various quasioscillations 
in the system. This can still better be seen by inspecting the connection of Aab((o) with 

(r) 

the G^g((o). We start from the observation that further complicated singularities may occur 
in the co-complex plane on the Riemann sheet of the analytical function Gab(g>) below the 
real axis where G^'2(co) is not an analytical function^. Then, in order to see what may hap- 
pen, let us assume that, in this region, at least approximately, a complex pole of the type 

(r) 

G)ab = G>AB — ^rAB (rAB > 0) may take place for G\g{(i)). In this simple situation, one can 
write 

,(r) ^ Cab 

(0 - (Oab + 'Tab 



GXilco) - „ „ ..^ (2.39) 



The advanced GF G)^g (co) has not a direct physical meaning and hence its possible singularities above the 
real axis will be not considered. 
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for real CO very close to the singularity value cOab- In Eq. (I2.39t . Cab is an inessential 
constant depending on the nature of the dynamical variables A and B. Besides, we suppose 
also that Aab((o) is a continuous real function of co so that we can use the exact relation 
(I2.26t . Then, from Eq. (I2.39t . we obtain the (approximate) expression 

TabRc Cab - (w - (OAB)Im Cab ,^ . 

Aab(0^2 ' . (2.40) 

((0-(0AB)2 + rj5 

If AAB(tt>) > 0, we should have necessarily Im Cab = and hence Aab((o) should have a 
Lorentzian shape. Unfortunately, according to the exact spectral decompositions (I2.32t and 
(I2.36t . Aab((o) is not a definite positive function in the whole range of real co. Hence, in 
principle, for each polar singularity for Gjj^(co) one must assume the quasi-Lorentzian or 
modified Lorentzian expression (I2.40t which reduces to a S-function singularity for co close 
to cOab and Fab —>■ 0. In any case, using the representation (I2.8t for 6(x), from Eq. (I2.39t 

(r) 

and the inverse Fourier transform for G^g(co), the time behavior of the retarded GF is given 
by 

Gi'2(T) = 6(1) ({A(x),B}) « e(x);?AB(x)^'-''"-^ (2.41) 

where 

j?ab(x) = ^ABe-^''\ (2.42) 

with SIab an inessential constant. Similarly, using the exact relation (12. 5t and the approxi- 
mate expression for ({A(x),B}) which can be extracted from Eq. (I2.41t . one easily finds: 

(A(x)B) - (A) {B) ~ CAB{^)e-'"'^'' = {CABe-^''^-''^''^ (2.43) 



Eq. (I2.43t is based on the boundary condition lim (A(x)B) = (A) (B) related to the physical 
feature that, for all real systems, there is always attenuation of the correlations on time. 
This is also supported by the rigorous Riemann-Lebesgue lemma in theory of analytical 
functions |23| as applied to the spectral representation (I2.7t where the spectral intensity 
Fab(co), and hence Aab(co), is a continuous function. 

(r) 

In conclusion, when Gjj^(co) has a pole below the real axis, with co close to cOab and 
Fab/ I COab I ^ 1> Aab(co) has a pronounced quasi-Lorentzian peak at co = COab and both 

(r) 

G\g{(i)) and (A(x)B) — (A) (B) are characterized by damped oscillations with frequency 
COab- Here Fab, which is related to the width of the peak, measures the damping of os- 
cillations. Of course, in this picture, Xab = ^ab defines the life-time of the oscillations in 
the system. It is also worth noting that, if Fab = (absence of damping), AAB(ro) has a 
5-function peak and undamped oscillations with frequency cOab occur in the system. 

The considerations made above for a single pole for G^^l{(i)) very close the real co-axis 
from below can be easily extended to several singularities of this type. In such a case, 
Aab((o) will result in a superposition of modified Lorentzian peaks whose widths measure 
the damping of the related oscillations. If the widths reduce to zero, Aab((o) becomes a 
superposition of 8-functions and this signals the occurrence of undamped oscillations in the 
system. 



142 



A. Cavallo, F. Cosenza, L. De Cesare 



In the next section, where we will consider the problem to search reliable approxima- 
tions within the classical GF formalism, we will see that the quasi-Lorentzian shape for 
A^b((o) is not mathematically adequate and another more appropriate representation, cap- 
turing the essential physics, will be necessary. 

3. Methods of Calculation 

All the previous general considerations suggest that, also in classical statistical mechanics, 
the formalism of the two-time GF's may be a useful tool to study the macroscopic prop- 
erties of classical many-body systems. The basic ingredients are Gas(co) and Aab(co) for 
two properly chosen dynamical variables A and B. Remarkably, if we are able to calculate 
Gab{c^) (Aab(co)), the exact relations established in the previous sections allow us to ob- 
tain Aab((o) (Gab((o)). In the present section we present the classical formulation of two 
general methods in strict analogy to the quantum counterpart, i.e. the classical EMM, for 
calculating the GF's, and the CSDM for a direct calculation of the SD. In principle, both 
the methods should give exactly the GF's and the related SD's and hence, in this sense, they 
are completely equivalent. Nevertheless, in practical calculations, previous experiences in 
quantum many-body theory E| |5l |2l IH IH QOl El El and, although reduced, in classical 
many-body theory ll2lfT3lfT4llT5ll22ll27l . suggest that the SDM has several advantages to 
make more reliable nonperturbative approximations in a systematic and controllable man- 
ner. 

3.1 The Classical Equations of Motion Method 

Differentiating Eq. (I2.2t with respect toz = t — t' yields 



Then, using the EM (12. 3t for dynamical variables and the obvious relation d%{x)/dx = 
5(x), Eq. (I3.lt becomes 



This is the basic EM for the GF ((A(x);B))y which, however, is not a closed differential 
equation. Indeed, in the right-hand side of Eq. (I3.2t a new higher order v-GF occurs 
involving Poisson brackets of a greater number of dynamical variables. Then, one needs 
to consider a new EM for the two-time v-GF (({A(t),//} ;B))^. The x-derivative of this 
function yields an additional equation, formally identical to Eq. (13. 2t with A (x) replaced 
by {A{x),H}, the right-hand side of which contains the new v-GF {{{{A{x),H} ,H} ;B))^. 
By iteration of this procedure, we obtain an infinite chain of coupled EM's for GF's of 




(3.1) 



^ ((A(x);B)), = 8(x) {{A,B}) + (({A(x),//} . 



(3.2) 
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increasing order which can be written in a compact form as 



^((x;,"A(T);B)), = 5(x)(WA,B}) + ((xr^A(x);B»^ , (m = 0,1,2,... 



). (3.3) 



Here Lh = iL = {...,H} and L'^A means L^A =A, l}j = {A,H}, Lj^A = {{A,H},H}, and 
so on. Notice that the chain of EM's (13.31) is formally the same for different types of GF's 
and hence one can eliminate the index v when the physical context is clear. 

In the practical calculations it is generally more convenient to work in the co-Fourier 
space. With ij^^dxe^^^ {df{x)/dx) = 0)f{o)), the chain of equations, in the co-representation, 
assumes the form: 



(0((i:;^'A(T);B)),,„ = /({x;^'A,B}) + /((i:-+U(T);B»^^„ , (m = 0, 1,2, ...). (3.4) 



To solve the chain of EM's in the form (13.31) or (I3.4t . we must add appropriate boundary 
conditions. For this, it is again more convenient to use the co-representation (I3.4t for which 
they can be obtained in the form of spectral representations or dispersion relations for the 
GF's (see Subsec. 2.1 and, in particular, Eq. (I2.18t ). in strict analogy with the quantum case. 
Anyway, although Eqs. (13.31) and (13.41) are exact, a free solution for interacting systems is, 
of course, impossible. In practical calculations one is forced to introduce decoupling proce- 
dures, and hence approximate methods, to reduce the infinite chain of coupled equations in 
a finite closed one which, although approximate, may be solved. However, in general, sys- 
tematic and controllable decouplings are not easy to find and one must check the reliability 
of a given approximation by comparing the results with experiments, simulations or other 
types of approaches. The CSDM, which will be the subject of the next subsection, seems 
more flexible in such a direction. 

3.2 The Classical Spectral Density Method 

Here, we present a general formulation of the CSDM fT2lfT4l for a systematic calculation 
of the CSD. For this aim, it is convenient to start from the expression (I2.11t of the SD in the 
T-representation. By successive derivatives of Aab(x) with respect to x and using the EM 
(12.31) . one has 




(3.5) 



(-/co)'"Aab(co) = / 



^T^''"^({X-A(T),B}). 



(3.6) 



Finally, integration over CO yields: 




(3.7) 
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The quantity on the left-hand side of Eq. (I3.7I) is called the m-moment of Kab{'^) and 
the relations (I3.7t constitute an infinite set of exact ME's or sum rules for the CSD. Notice 
that, for m = 0, the sum rule (I2.15I) is reproduced. 

Equations (I3.7t can be seen in a different way. Indeed, since the Poisson brackets 
and hence the ensemble averages involved in the right-hand side of relations ( 13. 7t can be 
evaluated, at least in principle, it is quite remarkable that the m-moments of the SD can be 
explicitly obtained without the a priori knowledge of the function A^b((o). This important 
result implies that the sequence of Eqs. (13.71) represents a typical moment problem. Its 
solution would yield the unknown SD and hence all the related quantities (GF's, correlation 
functions and other observables). Unfortunately, also this problem, in general, cannot be 
solved exactly and one must look for an approximate solution. Suggested by the exact 
classical spectral decomposition (I2.32t or (I2.35t . we seek for an approximation for Aab((o) 
as a finite sum of properly weighted 5-functions of the form (polar approximation) 

Aab((o) = 271 'kk{A,B)h i{(0-(Ok{A,B)) , (3.8) 

k=\ 

where ?i is a finite integer number. The unknown parameters ?i/:(A,B) and (0| (A,B), depend- 
ing on the nature of the dynamical variables A and B, have to be determined as a solution of 
the finite set of 2n (generally integral) equations obtained by inserting the expression (13. 8t 
in the first 2n ME's (I3.7t . This is the basic idea of the CSDM. Physically, the parameters 
(i)k{A,B) in Eq. (13. St play the role of effective eigenvalues of the Liouville operator and 
each of them, as a real (approximate) pole of the GF Gab((o) (see the spectral represen- 
tation (I2.20t or (I2.37t ). corresponds to a possible mode of undamped oscillations for the 
correlation function (A(x)B) (see Eq. (12. 38^ ). 

As discussed at the end of Sec. 2, there are physical situations where the damping 
of classical oscillations in the system under study may be important and hence the polar 
approximation (13. 8t is inadequate. In these cases, the basic idea of the SDM, related to the 
moment problem (I3.7t . remains still valid but it is necessary to chose a more appropriate 
functional structure for the SD which allows us to determine the modes of oscillations in 
the system and their damping or life-time. An extension of the SDM in this sense was first 
proposed by Nolting and Oles [28 J for Fermi systems and by Campana et al. [ 13 1 for Bose 
and classical systems whose SD's are not positive definite in the whole range of co. Here 
we will focus on the classical case. In Sec. 2, we have shown that information about the 
spectrum and the damping of the oscillations in the system can be obtained if the SD is 
assumed as a superposition of the quasi-Lorentzian peaks. The peaks would correspond to 
the frequencies of oscillations and the widths would measure their damping or life-time. 
However, we are now in a position to easily check that modified Lorentzian shapes for 
Aab((o) cannot be valid over the whole co range since all the SD moments of order m>2 
would diverge. This feature clearly contradicts Eq. (I3.7t . To assure the convergency of 
the SD moments at any order and to preserve the intrinsic physical character of the CSD's 
connected with the not everywhere-positive factor (O in the exact spectral decomposition 
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(I2.32t or (I2.35I) . one can assume the modified Gaussian ansatz for the Aab((o) IITbI 

exp 

,=1 V^r,(A,B) 

Notice that, with the functional representation ( I3.9t for Aab((o), the parameter ri^(A,B) 



Aab((o) = 271(0 £ ?i^(A,B)- 



-(m-a)t(Afi))^ 
r*(A,B) 



(3.9) 



has not to be identified directly with the parameter Y{A,B) in Eqs. (12.391 ) for Gab((o) and 
(12.401 for Aab((o) but, rather, with r^(A,B). Then, the width of the peak in co = (0<;(A,B) 

1 12 

is represented here by (A,B), the condition F(A,B)/ |a)(A,B)| ^ 1 must be replaced 
by Tk{A,B) / [(iik{A,By^ <^ 1 and the life-time of the classical excitations with frequency 

1 /2 ^ _ _ _ 

(Ok{A,B) is identified by 'lkiA,B) = F^' (A,B). The choose ( 13. 9t is only motivated by the 
fact that it takes direct contact with the notation used in literature fT3l l28l and simplifies 
the algebra in explicit calculations. Here the parameters (0/;(A,B) represent the oscillation 
frequencies and the new parameters Fj.(A,B), which describe the broadening of the 8-poles 
due to the finite life-times of the respective oscillations, give a measure of the damping 
effects. In any case, the representation (I3.9t assures that Aab((o) is a superposition of 
sharp quasi-8-f unction peaks and reduces to the polar one (13. 8t in the limit rk{A,B) 0, 
as expected. This functional structure allows us to investigate classical oscillations and 
damping effects systematically just as the 5-function ansatz (13. 8t does for case of undamped 
oscillations. Of course, with the ansatz ( 13. 9t . the parameters ^/({AjB), (s)k{A,B) and rif{A,B) 
have to be calculated solving the first 3n ME's (I3.7t consistently with the basic condition 
F,(A,B)/[(0,(A,B)]2«1. 

As in the EMM, also in the SDM the problem remains to close the truncated finite set 
of ME's arising from the polar ansatz (13. 8t or the modified Gaussian ansatz (13. 9t . In any 
case, the evaluation of the right-hand sides of Eqs. (13. 7t generally involves higher order 
SD's. Therefore, higher order moment problems should be considered, but the difficulty of 
calculations will increase considerably. So, in order to solve self-consistently the finite set 
of ME's, which arises from Eq. (13. 7t using the ansatz (13. 8t or (I3.9t . it is usually necessary 
to use some decoupling procedures and thus to introduce, in a systematic way, additional 
approximations in the CSDM which, however, will be automatically consistent with an 
increasing number of sum rules for the SD. 



4. The Classical J-dimensional Spin-5' Heisenberg Ferromagnet 
with Long-Range Interactions: A Many-body Approach 

4.1 An Introduction to the Model 

The classical J-dimensional spin-5 Heisenberg ferromagnet with LRI's decaying as r^P 
(p > d) with the distance r between spins, in the presence of an external magnetic field, is 
described by the Hamiltonian 

^ = -^ I JijSi-Sj-hY^S]. (4.1) 
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Here A'^ is the number of sites of a hypercubic lattice with unitary spacing, 
^Sj = {SpS-'pS^j); j = 1, ...,A'^| are the classical spins, h is the external magnetic field and 
the exchange interaction, in view of the thermodynamic limits as A/^ — > oo, is assumed to 
be Jjj = J/rfj, where rjj = |r,- — rj\ and J measures the strength of the spin-spin coupling. 
The extreme case p —^oo corresponds to the standard nearest-neighbor interaction while the 
mean field approximation is obtained when p = (replacing / by J/N). For this type of 
interaction, the thermodynamic limit A/^ — > oo is well defined only for p> d, while for p <d 
the ground state of the system has an infinite energy per spin asN and the conventional 
statistical mechanics cannot be directly applied. 

The classical spin model (I4.lt can be properly described by the set of 2N canonical 
variables {(pj,S^j} where cpy is the angle between the projection of the spin vector Sj in 
the {x — y)-p\ane and the x-axis. The Poisson bracket of two arbitrary classical dynamical 
variables ^(cp,^^) and B(cp,5^) is then defined by 



N 



dA dB dA dB 



For the following calculations, we find convenient to introduce the new spin variables Sj^ = 
S^- ± iS^j, so that 5^ = (S^j)^ + SjSJ . Then, if one defines the Fourier transforms of the spin 
vectors and the exchange interaction as 

^J = I^"'''^; ' J{k) = Ze-''<^^--^^'bij, (4.3) 
7=1 7=1 

where k denotes a wave vector in the J-dimensional Fourier space, the Hamiltonian (4.1) 
can be rewritten as (with j(k) = J{—k)) 

^ = -^Y.m -M^. (4.4) 

The sum in Eq. (14.41) is restricted to the first Brillouin zone (IBZ) of the lattice. Within this 
representation, the Poisson brackets for the spin Fourier components relevant for us are 

{st,Sl,} = ±iSl^^, , {s^,S,} = -2iSl^. (4.5) 

On the basis of the previous definitions, it is now easy to establish exact general formu- 
las for the internal energy and the free energy per spin from which all the thermodynamic 
properties of the spin model (14.11) can be obtained when the appropriate GF's or SD's have 
been determined using the many-body formalism presented before. For this aim, the Fourier 
representation ( 14.41 ) of the Hamiltonian in terms of the dynamical variables StS^^ and SiS^ ^ 

k —k k —k 

is particularly convenient since, in general, all the thermodynamic quantities can be easily 
expressed exactly in terrr 
we rewrite Eq. (14.21) as : 



expressed exactly in terms of the correlation functions ( S^^S ^ ) and ( 555'' ^ ) . To see this. 



H = Hq + Hi{J), (4.6) 
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with 
and 



Ha 



-hSL 



Hi{J) 



2N 



(4.7) 
(4.8) 



N 



where y(k) = j(k)/J = I e'*^ (^'-^^V''o- Then, the internal energy of the spin model is 
formally given by the ensemble average 



where 
and 



U = Uq{J) + Ui{J), 



Uo{J) = {Ho) = -h{Sf,), 



UjiJ) = {Hj)iJ) = - — l^y{k) 

1 

For the free energy F = — P^UnZ, we have 

Then, by integrating Eq. (I4.12t between and J with the initial condition F{J = 0) = Fq 

-p-Hn [Jdre-^"°],we get 



(4.9) 
(4.10) 
(4.11) 

(4.12) 



(4.13) 



From these general expressions, it is easy to obtain, for the internal energy u = U /N and 
the free energy f = F /N per spin, the desired exact relations (with N —^°°) 



(4.14) 



and 



k 







k -k/j' \ k -k/ji 



(4.15) 



Here, /o is the free energy per spin for a magnetic model without interactions and 
denotes an ensemble average as function of the interaction strength /. So, when the cor- 
relation functions (^StS^S^ and (^SiS^ are known, Eqs. (14.141 and ( 14.151 ) allow us to 
determine u an f and hence all the thermodynamic quantities using standard relations. 
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Also the transverse (^j^) and longitudinal (^j) correlation lengths can be obtained in 
terms of the previous exact relations. This can be performed by using the definitions 



It is worth nothing that, since the SDM allows to calculate (^Si^S and (^S^S^ solving 
integral ME's, all the thermodynamic quantities can be obtained without the explicit calcu- 
lation of the partition function. This is a relevant aspect of the SDM as a moment problem 
both in quantum and classical statistical mechanics. 

As a conclusion of this subsection, to fully appreciate the effectiveness of the CSDM, it 
may be useful to summarize the main known features about the model (14. It - (l4.4t obtained 
by means of rigorous methods, MC simulations and alternative microscopic techniques. For 
a future comparison, we include also some information for the corresponding quantum spin 
model which attracted recently much interest. 

The one and two dimensional long-range quantum spin- 1/2 Heisenberg ferromagnets 
(for brevity reasons, models of the type (I4.lt will be also named long-range spin models) 
in absence of an external magnetic field were investigated by Nakano and Takahashi using 
the so called modified spin-wave theory ll29l and the Schwinger-boson mean field approxi- 
mation l30l . Further information were derived for the cf-dimensional case by means of the 
EMM for the two-time GF's using the Tyablikov decoupling procedure [31] (in the next 
section we will present shortly the corresponding classical version for model (I4.lt ). Monte 
Carlo simulations for the two-dimensional quantum spin- 1/2 Heisenberg model have been 
performed for 2 < /:) < 6 ll32l . This scenario has been recently enriched by an extension 
||20| of the Mermin- Wagner theorem ll33l for the existence of ferromagnetic (FM) long- 
range order (LRO) at finite temperature in quantum Heisenberg and XY models in d{= 1 , 2) 
dimensions with r^P— and oscillating— exchange interactions. 

Classical long-range spin-5 Heisenberg FM models have attracted great attention, too. 
It has been rigorously proved the LRO exists in d(=l,2) dimension when d < p <2d \16 17 1 
and is destroyed at all finite temperatures for p>2d llT8l[T9 l. Similar results were obtained 
for the spherical model [341 and the present scenario of the critical properties is largely 
based on renormahzation group calculations for the classical n-vector model |35l|36l. Clas- 
sical long-range antiferromagnetic (AFM) models have been studied less extensively. The 
available rigorous results P16l ll7lfT8l suggest orientational disorder at all finite temperatures 
when p > 2d, but no theorem exists entailing existence or absence of LRO ford<p< 2d. 
Monte Carlo simulations have been also performed for both d(= 1, 2)-dimensional classical 
FM (for p = 2d 1211 ') and AFM (for p = ^ and p = I with d = I and d = 2, respectively 




(4.16) 



and 




(4.17) 
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f22\) Heisenberg long-range models. The results confirm that FM-LRO survives at finite 
temperature provided d < p <2d and allow to conjecture that no AFM-LRO exists at all 
finite temperatures for p> d. Spin-wave studies ll22l agree with the last conjecture but no 
definitive statement can be drawn at the present stage. 

Recent studies for model (14. U |24, 27 1 via the CSDM to lowest order approximation 
are in good agreement with previous available analytical and numerical investigations. Be- 
low we present the problem in a detailed and more systematic way with the inclusion of 
additional results showing the efficiency of the SDM also in classical statistical mechanics. 
In next subsection, as mentioned before, we present also the essential results for the model 
(I4.lt using the classical EMM and the Tyablikov-like approximation, although our main 
interest is devoted to the CSDM which will be used extensively in the remaining part of the 
article. In any case, these results will be very useful as a comparison and to underline the 
potentiality of the classical many-body formalism presented in Sees. 2 and 3. 

4.2 The Classical Equations of Motion Method and the Tyablikov-Like 
Approximation 

For the Hamiltonian model in the form (I4.4t . we introduce the GF (without distinction 
between the retarded and advanced one) 

G^,{t-t') = [(sl{t);S-_.^{t'))), (4.18) 

with the Fourier transform 

= ((5+(x);5>J)^ = j^JxG^^{x)e''^ . (4.19) 
From Eq. (I3.2t . we have for Gj((o) the EM 

coG,(a))=/({4,5:J) + /(({5^+(x),//};5>J)^. (4.20) 
Then, using the spin Poisson brackets (14. 5t . Eq. (I4.20t reduces to 
((0-/z)Gj((o) =2Nm + 

p 

(4.21) 

where m = (^S']^ is the magnetization per spin. This equation, as explained in Subsec. 3.1, 
is not closed and hence a decoupling procedure is necessary to determine G^((o). In strict 
analogy with the quantum Tyablikov decoupling |2|, we introduce the classical Tyablikov- 
like decoupling (or random phase like approximation) for which one can approximate 




(4.22) 
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Then, Eq. (14.2 It yields for Gj((o) the closed equation 

((0 - h)G7{(0) = 2Nm + m (j{0) - J{k)) Gv((0 



(4.23) 



This gives {{Ty) stands for Tyablikov): 



G 



CO) 



2Nm 

(Ty) 

k 



(4.24) 



where 



(Ty) 

k 



h + mjh{0)-y{k)) , 



(4.25) 



is the Tyablikov-like dispersion relation. From Eq. ( I4.24t and the 5-function representation 



lim — 
e^0+ 271/ 



1 



1 



x — ie x + ie 



5{x) 



(4.26) 



one immediately has (see Eq. (14.211) ') 

a('~>')((0) = / [Gj((o + /£) - G^((0 - /£)] = 2n{2Nm)5{(i) - 4^''^) 



(4.27) 



Notice that the expression (I4.27t for the spectral density a[^-*^^ (co) satisfies the sum rule (see 

Eq. dnu) 



2Nm. 



(4.28) 



Finally, from the spectral representation (I2.14t . we find the time-dependent correlation 
function /st{x)S^^) to be (hereafter we put Kb = I) 



s^{.)s:^ 



(Ty) 



,/a)Ap'((o) . 

271 CO 

(Ty) 



2NmT _i^y)^ 
-e * . 



CO 



(Ty) 



(4.29) 



(Ty) 



In the static case (x = 0), one has (^S^S -.j = 2NmT /ox. 

The previous relations contain the unknown mean value m = (S^) and hence, in order 
to determine the thermodynamic properties of our spin model, we must obtain a suitable 
expression for m. This constitutes a serious difficulty for classical spin systems because 
one cannot write, also for 5 = 1 /2, a classical counterpart of the quantum kinematic rule 
for the z-component of the spin. However, working in the low-temperature regime where 
the angular momentum is nearly saturated {S^j ~ S), from the identity 5^ = {S'j^ + SjSJ 



we have 
in the form 



SjSj 



/S^ and the magnetization per spin m can be approximately expressed 
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Then, from Eq. (I4.29t . one should have: 



m ~ 5 < 



Tm ^ 

k »j 



{Ty) 



(4.31) 



Equations (I4.25t and (14.3 It constitute a closed self-consistent system for (co^ , m) which 
can be solved to obtain the thermodynamic properties under near saturation condition. 

A general and reliable expression for m in terms of the dispersion relation, valid for 
arbitrary temperature and including the near saturation expression ( I4.31t . was suggested in 

Refs. [14 24 J . However, since our main purpose is to focus on the CSDM predictions, 

we postpone the derivation of such general expression to the next subsection. Then, work- 
ing to the lowest order in the CSDM, we will obtain a self-consistent system of equations 
which differs from that in the Tyablikov-like approximation only for a new expression of 
(0^. Hence, once these more general equations are studied, one can similarly discuss the 

classical Tyablikov ones simply replacing co^ with co^^'^^- 

4.3 The Classical Spectral Density Method and the Moment Equations 

Let us introduce the co-dependent CSD 

A,(a)) = /({^ (x),^:,})^ = //;;^xe-({4(x),5:J) . (4.32) 



The system of ME's for A^((o), with A = and B = S - in Eq. (I3.7t . is given by 

/7^»"'^?(») = -^'""'({^^4'^:J) ' ('« = 0,l,2,...). (4.33) 

Focusing on the first three ME's for next considerations, the use of the basic Poisson brack- 
ets (l43l yields 



N . 
k 



(4.35) 
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£ \j(ki)-j(ki+k)\ \j(ki-k2)-Jik2) 



{ 



k^—ki —ki k[ 



ko—ki —k—k2 k-\-ki 



ki—ki k[ —k2 } ~^ 



2h- 



(4.36) 



where we have made frequent use of the standard properties of the Poisson brackets. 

Notice that, in view of the general Eq. \2Ab\ . in the previous ME's only the correlation 



function {StS 

k —k 



can be exactly expressed in terms of Aj((o) as 



' <i(oAj((o) 
271 Poo 



(4.37) 



As explained in Subsec. 3.2, in order to give an approximate solution to the "moment 
problem" (14.331) (or (I4.34t - (l4.36t ). according to the CSDM one must choose an appropriate 
functional structure for the SD. The aim is to truncate the system of moments (I4.33t at 
different levels and to determine A^((o) as a solution of a finite set of equations. In the 
next sections we will consider some of the numerous possibilities suggested in Subsec. 3.2. 
Although the chooses are very simple in order to avoid prohibitive calculations, they are 
able to capture the essential physics of the spin model under study. In any case, one must 
remember that finding Aj((o) as a solution of a finite set of ME's means to have a SD which 
satisfies two or more sum rules. 



4.3.1 One 5-function Polar Ansatz 



As suggested by the exact spectral decomposition (I2.35t . according to the spirit of the 
CSDM, we try to determine Aj((o), to the lowest order, in the form llT3ll24ll27l 



A7((0) =27tXr5((0-(Or), 



(4.38) 



involving two unknown parameters and co^. Then, to determine these parameters only 
the first two ME's (I4.34t and (I4.35t are necessary. Taking into account the polar ansatz 

(4.39) 



(I4.38t . these reduce to the integral equations: 



2Nm 



h'^k 



1 



2Nmh + -J^ J{k')-J{k-k'] 
k' 



StS:^ +2(S^S^ 



(4.40) 
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To close this system, we must express the involved unknown quantities in terms of A7((o). 



The transverse correlation function (^StS -^j, which enters the right-hand side of Eq. (I4.40t . 

can easily expressed in terms of the SD and hence in terms of the parameters and coj. 
Indeed, Eq. (12. 16t and the ansatz (14.381) yield 



Sp-^ = ^. (4.41) 



Unfortunately, there is not a simple way to express the longitudinal correlation function 
{s^S^- in terms of Aj((o). Hence, to close the truncated system of ME's (I4.34t - (l4.35t . 

one should introduce the longitudinal SD a1^^((&) = i(^^S^^{x),S^ and then formulate 
and solve another moment problem according to the basic idea of the CSDM. However, 
this problem would be coupled to that for the transverse SD A^((o) and the difficulties 

would be sensibly amplified also if a one 8-function structure for both Aj((o) and A^^^((o) 
is assumed. The simplest way to overcome this difficulty, frequently used in literature 
El|5l|2l|8l|9l[ini[IIl[nili3lIIl[l5|,isto resort, as a first level, to the decouphng procedure 
( ^1^^ l) ~ ^^l) l) —^^^^^ which corresponds to neglect the correlations between 
the Fourier transforms of the longitudinal spin components. With this approximation, Eqs. 
(I4.39t - (l4.41t immediately yield the self-consistent equation for the frequency dispersion 
relation 

^, = h + m (7(0) -m) + ll^ J{k')- J{k-k') 



It is worth noting that the first two terms in the right-hand side of Eq. (I4.42t constitutes just 
the expression for the dispersion relation obtained in Subsec. 4.2 using the Tyablikov-like 
decoupling within the EMM. Of course, the previous decoupling and the Eq. (I4.42t are 
appropriate to describe only thermodynamic regimes with a finite magnetization as under 
near saturation conditions. In particular, one cannot use the basic Eq. (I4.42t to explore near- 
zero magnetization domains in the phase diagram as the critical region and the paramagnetic 
phase in zero external magnetic field. To overcome this difficulty, one is forced to find 
a more appropriate decoupling procedure for the longitudinal correlation function which 
allows us to obtain self-consistent ME's appropriate for describing regimes when m — > 
and preserves also the simplicity of the one 5-function ansatz for the transverse SD. A 
possible and successful solution to this crucial question was suggested many years ago for 
Heisenberg spin models with short-range interactions (SRI's) |14|. It was shown that a 
suitable decoupling procedure when the magnetization approaches to zero (see also Refs. 
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141 1371 for the quantum counterpart) consists in writing in the ME (I4.35t 
1 



1 



£ Jik')-Jik-k') 



k' 




(4.43) 



This equation only involves Aj (co) (see Eqs. (I2.16t and (14.41^ ) and hence, inserting in the 
ME ( I4.35t one finds, for the dispersion relation, the new expression 



oyr=h + m(j{0)-J{k)) + 



T yJ{k')-J{k-k') 



^k' 



(4.44) 



Notice that the effect of the decoupling ( 14.431 ) corresponds essentially to perform in Eq. 
(l4!42l the transformation !(...) m^/S^Z{-)- Of course, alternative decoupling proce- 

k' k' 

dures can be conjectured for taking into account the effect of the longitudinal spin corre- 
lations working in terms of the single Aj((jo). An example, appropriate to near saturation 
regimes, will be presented in Subsec. 4.5.2. 



4.3.2 Magnetization per Spin in Terms of the Transverse Spectral Density 



To close the system of ME's (I434l i- (l436t for Aj((o) (and also the EM m2\\ for Gj((o)), we 



must find an explicit expression for m = {S'jj in terms of the SD. A reliable and successful 
expression of m, valid for any T and h, can be obtained by means of the following procedure 
within the spirit of the SDM fl4ll24l . Let us introduce the higher-order SD: 



kjc' 



(4.45) 



where the summand is constructed by A^((o) = i(^^Si{x),S with the change 5^ (x) 



5t 7,(x) and 5 



, ^ S -.S'' The SD (I4.45t is associated to the higher-order GF 

— K — K — r 



G('-)((0)=£ 

kl' 



e(x)({4^,(x),5:,5i,} 



(4.46) 



We now perform in Eq. (I4.45t the decoupling procedure 



{'Uk'^'^^'Zk'U'})^ - {'-k') (1 ({?+^'(^)'^: J) ' (4.47) 
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where a is an unknown parameter to be properly determined (notice that, with a = 0, Eq. 
(I4.47t reduces to a Tyblikov-hke decoupling). It is easy to see that this allows us to reduce 
n((o) to the simplest form 



a((0)=m(l+fl)-£A 



(0 



(4.48) 



in terms of the original Aj((o). The parameter a can be now determined by imposing a 
sum rule for n((o), and precisely, by requiring that the zeroth-ME for the exact n((o) is 
preserved in the decoupling procedure. Taking into account the general Eq. (12. 15t and the 
identity {S^jf = 
for a((o) in Eq. (l4!45l . 



-a((o) 



SjSj , a straightforward use of the Poisson bracket properties yields 



2N^S^ -3N(SpJ 



(4.49) 



On the other hand, the zeroth-ME for the reduced form ( 14.481 ) for n((o) is given by (see Eq. 

^ — a((0)«m(l+a)-X: y_ — Aj((0)=2A^m2(l+a). (4.50) 



(14341 ) 



Then, by imposing that Eqs. ( I4.49t and ( I4.50t coincide, we obtain 

m\l+a)=S^-^(^S+SJ 



(4.51) 



From this equation, it follows that in the high-temperature regime, when the isotropy 

3S^/2 is consistently reproduced regardless of the value of a. On the 



condition (SjSj 



other hand, in the regimes where the spins S^j are nearly saturated (S'j ~ 5 or SjSj /S^ <^ 
1), the magnetization per spin can be approximately expressed in the form ( 14.301 . Then, 
according to Eq. ( I4.5U . the parameter a must go to zero in an appropriate way so that Eq. 
(I4.30t is recovered. Bearing this in mind, a direct comparison between Eqs. (I4.30t and 
( 14.5 It provides for a the expression 



2S2 



(4.52) 



Inserting it in Eq. (14.511 ) leads to the required formula for the magnetization per spin m 
suitable for our aims (with h> 0) 



m 



1-- 
2 



3 \SJSJ 



52 



'-2 



1 



52 



(4.53) 
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where 



SjSJ 



1 



(4.54) 



One can immediately check that, under near saturation condition, Eq. ( 14.531 1 reproduces 
.30t . In view of Eq. (14.541 and the exact expression of 



the relation 



in terms of 



A^((o) (see Eq. (I4.37t ). Eq. (I4.53t for m allows us to take immediate contact with the SDM 
and, in particular, with the one 5-function for A^((o) introduced in the previous subsection. 

Postponing this problem to the next subsection, we wish to outline here how the expres- 
sion (I4.53I) . which is valid for arbitrary temperature and magnetic field, can be used in the 
EMM for Gj((o) within the Tyablikov-like approximation (see Subsec. 4.2). In this case 

the GF has a real pole at d^^ '' = h + mJ{y{0) — y(k) which corresponds to a one 6-function 
for the SD. This feature and the general expression ( I4.53t for m allows us to express m as 
a function of d^^'^ yielding to a system of two self-consistent equations. When this system 
is solved, all the thermodynamic properties can be obtained. 

It is worth emphasizing that the one 5-pole ansatz for A^((o) within the CSDM, involv- 
ing two ME's, is not equivalent to the classical Tyablikov decoupling. This becomes clear 
by comparing the two different expressions (I4.25t and (I4.42t for the dispersion relation, 
although one assumes the same expression (I4.53t (or (I4.28t ) for m. The main reason is 
4^-^V(o) (and hence col"^-^'''' 

k k 

Aj (co) (and hence coj) is obtained to satisfy both the sum rules (I4.34t and (I4.35t 



that, while AV (co) (and hence coi ) satisfies only the first of the sum rules (I4.34t - (l4.35t . 



In this 

sense, one can claim that the one 5-function ME's solution for A^(co) in the CSDM is bet- 
ter than the corresponding one derived by means of the Tyablikov-like decoupling in the 
classical EMM for Gj(co) and constitutes a one step beyond such an approximation. 



4.3.3 Moment Equations and Dispersion Relation for one 6-pole Ansatz 



We now come back to the Subsec. 4.3.1 and add, to the Eq. (I4.42t or (I4.44t for the dis- 
persion relation, the corresponding equation for m in terms of the parameter and coj in 
A^(co) = ?ij5(co — coj). From Eq. (I4.53t . we easily have 



1-3— y- 

S2 ^ Tniy I 

k k 



which, in the near saturation regime, reduces to 



(4.55) 



m '. 



Tm • - 1 



(4.56) 



Eq. (I4.42t (for m ^ 0) or (I4.44t (for m ^ 0) and Eq. (14.551) represent the closed system 
of ME's to be solved self-consistently. For this purpose it is convenient to introduce the 
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dimensionless variables 



m 

^ = ^ 



T 



JS 



(4.57) 



Then, taking the thermodynamic limit A/^ —> 00 with J^g^d'^k{...) / (27i)^, our self- 



constistent equations can be written as 



l-3ra 



IBZ {2tiY (Or 



l-3Ta 



d'^k 1 

isz (271)'' COj 



with 



m 



1+^ 



l + fa 



d'^k' Q.^P\k -k') - Q.'-P\k') 
iBZ (InY Q.(p){k)Gy^, 

d^'k' a^p\k-k')-a^p\k') 



IBZ (27l) 



and 



a(p\k)=y{0)-y(k)=l^ 



1 — cos(^-r) 



a/0 



0, 



(4.58) 
(4.59) 



(4.60) 



(4.61) 



The system of self-consistent integral Eqs. ( I4.58t - (l4.60t for the unknown parameters 
G)^ (the dimensionless or reduced excitation dispersion relation) and a (the dimensionless or 
reduced magnetization per spin) is very hard to solve. For obtaining explicit results one is 
forced to consider asymptotic thermodynamic regimes or resort to numerical calculations. 



4.4 Thermodynamic and Critical Properties within the one 6-pole Ansatz 
for the Transverse Spectral Density 

4.4.1 Low-Temperature Properties and Long-Range Order 

Our primary purpose is to solve the self-consistent system of Eqs. ( I4.58t - (l4.60t and then, 
using the general relations in terms of A^((o), to determine the relevant thermodynamic 
properties of the classical Heisenberg FM model ( 14. U . Since any attempt to obtain explicit 
analytical solutions for all the allowed values of T and h is hopeless, we first examine the 
possibility to have analytical results in the low-temperature regime. We can expand Eqs. 
(I4.58t and (I4.59t in power series of f for allowed values of h which prevent the occurrence 
of divergences. Focusing on the first order expressions in the reduced temperature, which 
are sufficient to capture the relevant features of the low-temperature physics of the model 
(I4.lt . we have: 

cbj ^ /j + a^P^ ik)-T^^F {h)D.^P^ (k) +J 'f^ (h,k)] + 0{T^) , (4.62) 
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and 



where 



and 



OP^l-T:F}''\h) + 0{f 



d''k 



iBz i2TiYh+aip)(k) 

d'^k' D.^P\k')-D.'^P\k- 



-k') 



IBZ (27l) 



(4.63) 



(4.64) 



(4.65) 



h + Q.(p){k') 

At first, we assume h so that no convergency problems occur for the ^-integrals in 
Eqs. ( I4.64t and ( I4.65t . Unfortunately, since one cannot obtain an explicit expression of 
the function il^'^' (k) in terms of elementary functions for k in the whole IBZ for arbitrary 
values of the interaction parameter p > d, the integrations in Eqs. ( 14.641 and ( I4.65t cannot 
be performed explicitly. Nevertheless, for a sufficiently low external magnetic field, the 
dominant contribution to the integrals arises from the low wave-vector excitations. Hence, 
one can obtain an explicit estimate of (O^ and a assuming the dominant behavior of D.^p^ (k) 

in the IBZ as k- 



0, provided that the /z-dependent coefficient of T in Eqs. (I4.62t and 
(I4.63t remain finite. Bearing this in mind, one can show I20ll2ni22ll30l that, for p > J, we 
have for Q.^p^ (k) the low-k expansions 

^ AdkP-'' + Bdk^ + 0{k'^) 



Q.^P\k) ~ 



p^d + 2 



(4.66) 



p = d + 2 



Here, the coefficients Aj, B4 and Q depend in a cumbersome way on the dimensionality 
d, the interaction exponent p and the wave-vector cutoff A related to the geometrical def- 
inition of IBZ. Their explicit general expressions are inessential for our purposes and will 
be omitted. However, for case d = 2, they will be explicitly presented in Subsec. 4.4.2 
where our analytical predictions near criticality are compared with those obtained by MC 
simulations. 

The most relevant analytical and numerical results for case d = \ will be considered 
explicitly after a discussion of the low-temperature properties for general d with p> d. 

Taking into account the asymptotic behavior (I4.66I) . the integral !F^''\h) can be easily 
estimated and we get, to leading order in f , the following near saturation representations 
for the reduced magnetization per spin 

/ p AdAP-^' 



a(r,/i)~ 1-r - < 



h 2F1I1 



Kd ^dkk 
Jo 



d-i 



p — d' 



h + CdkHn 



h 



h 2F1 



d d + 2 



h j 



d <p<d + 2 



■-d + 2 



p>d + 2 



(4.67) 
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where 2F]_{a,b,c;z) is the hypergeometric function, = /T{d/2) and r(z) is the 

gamma function. 

An explicit estimate of (O^ for small k can be obtained by means of an analogous but 

rather complicated and tedious calculation of ;F2^\h,k). However, since an explicit esti- 
mate of (0^ for arbitrary p and d (with p> d)is inessential for next developments, we avoid 
to present the related cumbersome result. Information about the reduced excitation disper- 
sion relation (Oj in the whole IBZ will be given below for (i = 1 varying the interaction 
parameter p > I. 

The low-temperature susceptibility % = {dm/dh)j can be now easily obtained from Eq. 

%/J we have 



(I4.67t . Indeed, for the reduced susceptibility % 



da 

dh 



h 2^1 2 



p — d p — d 



Kd / dkk 



h + CdkHa 



--2 ( d d + 2 Dd. 



d < p <d + 2 



p = d + 2 



p >d + 2 



(4.68) 



Of course, the above expressions have a physical meaning only for values of f and h such 
that the near saturation condition (a ~ 1) is assured and hence also for /j — > if the long- 
range order (LRO) occurs. In particular, if we write in general a ~ 1 — a(/j)f, this condition 
implies necessary a(Ji)T <^ 1 and hence Ta(h) becomes the natural expansion parameter in 
the problem. As already mentioned before, the integrals (I4.64t and (I4.65t . and hence the 
functions in Eqs. 



(I4.67t - (l4.68t . could diverge as h 
When this does not happen, Eq 



for particular values of p and d. 



I4.67t should yield a spontaneous magnetization at a finite 
temperature signaling that our spin model exhibits LRO. From the low-^ behaviors (I4.66t it 
is easy to see that the integral (h) converges in the limit /j — > only foT d < p <2d with 
d <2 and for p > li with d > 2. This means that for these values of p and d a spontaneous 
magnetization per spin mo{T) = ao{T)S, and hence LRO, exists at small but finite T with 



ao(r)~i-f:r/^)(o) 



where J^^^^ (0) is a finite quantity given by Eq. (I4.64t for h 



(4.69) 

0. Of course, an estimate 
For 



of this constant (depending on p and d) can be immediately derived from Eq. (14.67 
p >2d with d <2, the integral ^j'^'' (h) diverges as h so that no finite solution for a 
exists at f 7^ 0. Then, for these values of p and d, Eq. (I4.63t is satisfied only if a = 
and hence no LRO occurs at finite temperature. In Fig. 1 we present the (p — <i)-plane 
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version of the phase diagram of our FM model displaying the scenario discussed above. 
Here we also show the regions where, as we shall see in the next Subsec. 4.4.2: i) the 
system exhibits a critical behavior like for the Heisenberg model with SRFs and ii) the 
LRI's become effective. 




Figure 1 : Regions of the (p — J)-plane where long-range order (LRO) takes place for a d- 
dimensional spin-S' Heisenberg FM model with interactions decaying as r^^. The dashed 
line separates the domains where long-range interaction (LRI)- and short-range interaction 
(SRI)- regimes occur. The dark region corresponds to nonextensive thermodynamics not 
considered through this paper. 



We can conclude that our low-temperature results suggest that a transition to a FM phase 
at a finite temperature occurs in the regions of the (p — (jf)-plane where LRO takes place. 
In the remaining domains a different scenario happens with absence of a phase transition. 
These predictions agree with the recent extensions 1 20 1 of the well known Mermin theorem 
l33l to spin models with LRI's of the type here considered. 

Other low-temperature thermodynamic properties follow from the general exact ex- 
pressions (14.141) and (14. 15t . Within our approximations, the correlation functions in these 
equations can be expressed in terms of the transverse SD and hence the reduced internal 
energy m = u/JS^ and free energy / = f /JS^ per spin assume the forms 



-Ir^-Tof ^fl-W, (4.70) 

JiBz {2%Y coj 2 
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and 



f = fo-f 



iBz {2ny 



y{k) 







^Y(O) f'dfa^f). (4.71) 

2 JO 



In particular, in the low-temperature limit, we find 



u^T-h- -y{0) + 0{T^g{h)), 



(4.72) 



where the explicit expression of the term 0{T g{h)) can be obtained from Eqs. (I4.62t . 



(I4.63t and (I4.70t in a straightforward way. A similar low-J expression can be easily ob- 
tained for /. From the result (I4.72t . the reduced specific heat at constant magnetic field 

Q = (du/dT)jj = Ch/S is given by 



Cj,c^\ + 0{T^g{h)), 



(4.73) 



as expected for a classical spin model P38V We now focus on the case d = I with p > I 
for which explicit analytical results can be obtained 1 24 1 allowing a comparison with recent 
MC simulations and transfer-matrix predictions. In this case, in the thermodynamic limit 
N ^ o°, we use the transformation ~^ lo^ '^^{■■■)/'2^ = lo '^^{■■■)/'^ (due to the 



symmetry of the problem) and, for small k 11291 1301 131 



r 1 



1 — cos(^?i) 



'Ml 



1 

I 2 



1<;?<3 
p = 3 



(4.74) 



!;ip-2)k^ , p>3, 



where r\{p) = kT ' (p) / sm[n{p — 2)] and ^(z) is the Riemann's zeta function. In particu- 
lar, for = 2, we get the exact results ll39l 



a(2)m = -)t--, 0<k<2n 



(4.75) 



With the low-A; expansions (I4.74t . for the reduced magnetization per spin a, we have the 
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low-temperature representation (see Eq. (I4.67t ) 

'■''dk 1 



a{T,h) 



l-T 



71 h + a^.P){k) 

1 p ■KP^^r[{p) 



\-T I 



h 9 F, 1 



2-ri 
71 



>-l>-l' 2 h 



2 \k 



, l<p<3 
, p = 3 



7i^(/7-2)V 2h 



where the hypergeometric function, when h ^0 behaves as 

1 p nP^^r\{p)' 



2h 



(4.76) 



h 2^1 I 1 



p-\ p-\ 
' ri(p)7iP-2 



2h 



2(2 -p) 



71^ 



In 



1 + 



71"" 
2/^ 



1 



1 



2h 



, 1 <p <2 
, p = 2 

, 2<;7<3. 



(4.77) 



According to Eqs. (I4.76t - (l4.77t . if we write, as before, a = 1 — Ta{h), the nearly saturation 
condition imposes that Ta{h) <C 1. For instance, with p>3, this condition is expressed by 
f I Vh <^ 1 and hence the natural expansion parameter is f/ Vh. Of course, more accurate 
estimates of the integral in Eq. (I4.76t can be performed including higher-order terms in the 
low-fc expansions (I4.74t . A check of the reliability of the estimates (I4.76t can be obtained 
by calculating the reduced magnetization per spin to the first order in T with p = 2 when 
n'^' (k) has the exact expression (14.751) . In this case, for /j > 0, we have 



a{f,h) ~ 1 



71 



■An 



(4.78) 



For small h, this reduces to 



a(f,/j) ~ 1-^f In 



71^ 



1 + - 



2h 



(4.79) 
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which is just the result obtained from the hypergeometric function representation in Eq. 

71 assuming the low-k expression D.^^\k) (see Eqs. 

14.781) (or ( I4.79t ) shows that the near saturation condition for — > 



K7B with p = 2 and r\{2) 



1477411 and (147771 ). Eq 

is satisfied only if rin(l//j) ^ 1. It is worth noting that, consistently with the previous 
results for general d (see Fig. 1), from Eq. (I4.76t . provided that, for small values of h (for 
finite h no problem occurs), the f-expansion preserves its physical meaning, the following 
low-r features arise: 

(i) for 1 < 7? < 2, LRO occurs with a spontaneous magnetization 



OoiT) = 1 



-T + 0(T' 



(4.80) 



(2-7^)11^)^ 

(ii) for /? > 2 and h 0, the coefficients in the f expansions in Eqs. (I4.76t - (l4.79t 
diverge and no LRO takes place at finite temperature. 

From the previous equations one can immediately obtain also the reduced susceptibility. 
In particular, for the simplest case p = 2 we get, exactly 



271 



f + h 



71 

h /-ni 



1 



:ln 



4 ^" 2 



(4.81) 



Besides, the internal and free energies per spin, and hence other thermodynamic quantities, 
can be obtained from Eqs. (I4.70t and (14.711) . In particular, for reduced internal energy per 
spin, we find (with y(0) = ^{p)) 



T-h-h{p) + 0{T^g{h)) 



(4.82) 



and hence C/- = 1 + 0{f^g{h)) as for a generic d. 

Consistently with the previous results and exact theorems for the J-dimensional spin 
model, our analytical near saturation calculations for the classical long-range Heisenberg 
FM chain with \ < p <2 predict, in a transparent way, a transition to a FM phase at finite 
temperature (see Fig. 1). For p>2,Si thermodynamic scenario at finite temperature, similar 
to that for SRI's 1 14 1, takes place with absence of a phase transition. 

Additional information about the excitation dispersion relation and thermodynamic 
properties for the FM chain, within nearly saturated regimes beyond the linear expansion in 
Taih) <^ 1, can be obtained solving numerically the full set of self-consistent Eqs. (I4.58t - 
(I4.60t with respect to (O/t and a for J = 1 |24 1. We present here the most relevant results in- 
cluding comparisons with predictions obtained by means of different (and in a some sense, 
exact) methods. As we will see, all the numerical results confirm the previous scenario 
based on the estimates involving low-temperature expansions (for appropriate values of h) 
and the dominant contributions to the dispersion relation as ^ ^ 0. In Fig. 2, the reduced 
dispersion relation (bt is plotted as a function of k in the interval 0<^<7iforJ'=0 (Fig. 
2(a)) and T = 0.5 (Fig. 2(b)) a.th = 0.5 with selected values of the interaction exponent p. 
The reduced magnetization a and susceptibility % in terms of T at a fixed value of h with 
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some values of p are shown in Fig. 3(a,b) where, as a comparison, transfer-matrix results 
for the Heisenberg FM chain with SRFs {p = oo) 1 14 40 1 are also presented. The plots of 
the reduced magnetization per spin a as a function of h for some values of f with p = 1 .5 in 
Fig. 4(a) and for different values of p with f = 0. 1 in Fig. 4(b) are particularly meaningful. 




1.00 



Figure 2: Dispersion relation cO/t of an Heisenberg FM chain as a function of k in the interval 
[0,7l] for reduced magnetic field h = 0.5 and temperature T = (a) and T = 0.5 (b) with 
different values of the decaying parameter p. 




Figure 3: Temperature dependence of the reduced magnetization a (a) and susceptibility % 
(b) of an Heisenberg FM chain for h = 0.5 and different values of p. The dots represent the 
transfer-matrix (TM) results for SRFs (/? = oo). 

They show clearly our previous analytical finding that, for 1 < p < 2, a FM phase occurs 
at finite temperature. Of course, at this stage, our ME's do not allow to study regimes with 
a — > as /j — > and hence to find the transition point to the FM phase where \ < p <2. 
Information about this transition and the related critical properties will be given in the next 
subsection. The reduced internal energy u and the specific heat Q, evaluated from the 
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Figure 4: Reduced magnetization a of an Heisenberg FM chain as a function of h for (a) 
p = \ .5 and different values of T and (b) T = 0.1 and different values of p. 



solutions of the {d = l)-ME's (I4.58t - (l4.60t . are also plotted in Figs. 5(a) and 5(b) as 
functions of T for different values of p. In both these figures, as a further comparison, we 
have also reported the transfer-matrix results for the short-range chain fT4l l40ll and MC 
simulations recently obtained 1 2 1 1 for a classical long-range FM Heisenberg chain with 
p = 2. Finally, in Fig.6, we present the behavior of the transverse correlation length as 
a function of T for h = 0.5 and different values of p. Notice that a substantial difference 
clearly occurs for 1 < < 2 and p >2, respectively. Again, in the same figure, we have 
plotted also transfer-matrix results for a FM chain with SRI's 1141 l40l . 



4.4.2 Thermodynamic Regimes witli Near-zero Magnetization 



As discussed in Subsec. 4.3.1, Eq. (14.421) for the dispersion relation, or the first expression 
for R(k) in Eq (14.601) . is inadequate to obtain suitable results for thermodynamics regimes 
with near-zero magnetization (a —>■ 0). However, this case can be successfully explored 
using the self-consistent ME's (I4.58I) - (I4.59I) with the simple modification for R{k) given in 
Eq. (I4.60t when a ^ 0. Below, we show indeed that the modified ME's can be properly 
used for an estimate of the main critical properties of our classical spin model, when it starts 
to exhibit LRO, and of the low-temperature paramagnetic susceptibility in the remaining 
regions of the (p — cf)-plane when no LRO exists. 



A. Transition temperature and critical properties In the limit /j — > with a > 0, the 
system of Eqs. (I4.58t -( l4.59t . with R(k) appropriate to include the case a ^ 0, assumes the 
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Figure 5: Temperature dependence of the reduced (a) internal energy u and (b) specific heat 
of an Heisenberg FM chain for h = 0.5 and different values of p. The stars represent 
Monte Carlo (MC) results for p = 2 and the dots TM predictions for SRFs {p = oo). 



form EH 



m 



T 



1 + 



1 



a^p)(k) JiBz (InY 
3TQiT) 



Q.(p){k')R{k') 



\-TQ{T) 



where 



QiT) 



1 



(4.83) 
(4.84) 

(4.85) 
(4.86) 



'lBz{2Ky a{p)(]c)R(k) ' 

From Eq. (I4.85t we can immediately obtain an equation which determines the critical 
temperature Tc (when it exists). Indeed, the reduced critical temperature Tc (with Tc = 
JS^Tc) can be obtained setting o{T c) = 0. This yields the equation 



i-3r,e(r,) =0. 



(4.87) 



It is not simple to solve exactly this equation and one must resort again to numerical cal- 
culations. However, it is instructive to get an explicit estimate of Tc by calculating the 
integrals in Eqs. (I4.83t -( l4.86t taking for Q.^^^ (k) the dominant contribution as ^ ^ 0. Then, 
solving by iteration our self-consistent ME's, to first level of iteration we have 



and 



(0) 



(4.88) 



(4.89) 
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100 




Figure 6: Temperature dependence of the transverse correlation length of an Heisenberg 
FM chain for h = 0.5 and different values of p (with logarithmic scales). The dots represent 
TM results for SRFs (p = oo). 



where ir/^'(0) anda(r) are given by Eqs 
Q{T) explicitly known, Eq 



I4.64t (for h = 0) and (I4.85I) . respectively. With 
BTHTt yields for T,: = Tc{p,d) 



1 



(4.90) 



From this simple expression, one immediately sees that a critical temperature exists only 
when !fI''\o) is finite and hence in the regions of the {p — d)- plane where LRO occurs. 
Here we explicitly consider the one-dimensional case for 1 < < 2 and the two-dimensional 
one for 2 < 7? < 4 for which a comparison with recent alternative analytical and MC results 
can be performed. 

For (i = 1 and I <p <2 one finds ll24l 



Tcid = \) 



(4.91) 



In particular, for p , Eq. ( 14.911 ) yields Tc = JS^T c ~ %^JS^{2 — p) /4 to be compared 
with Ti^ J S{S+ 1)(2 — p)/3, obtained in Ref. fJT\ for the quantum case by means of 
the EMM, and with 75^(2 — p) derived by Kosterlitz |36 1 for the classical chain using 
a perturbative renormalization group approach. 

For the case d = 2 with 2 < p < 4, the coefficients of interest in the low-^ expansions 
(14:661 for (k) are given by 1 3 1 ] 



A2 = 2^-Pn^r-^{p)/sm[Ti{p - 2)/2] 



: 2P-% 



p-2 



2-p 1 



2-p 3 
2 '4 



,2<p<4, 



(4.92) 
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where ^(z,a) = Lr=o(" + '^) ~ '^^e generalized Riemann zeta function and ^(z) = ^(z,0) 



is the ordinary Riemann zeta function. Then, in Eq. (I4.90t . we have 

^1 ^^i^2"4^2T^'^"[''^^-^)/^] ' ^^-^^^ 



and hence 



Tid-D- (4-^)^'^^'' (494) 
ic[.d-l}- sin [T,{p - 2)/2] ■ ^^-^^^ 



In Fig. 7 the critical temperature of {S = ^)-models, for d = \ (Fig. 7(a)) and d = 2 (Fig. 
7(b)), is plotted as a function of p and compared with the corresponding results recently 
obtained, for the long-range quantum Heisenberg model by means of the EMM, within the 
Tyablikov decoupling, for both cases d = \ and d = 2fi\\, the modified spin-wave theory 
(291 and MC simulations |32|, for <i = 2. Our results appear to be consistent with those 
obtained for the quantum counterpart if we keep in mind that the apparent discrepancy is 
related to the known feature that a quantum spin-S' model can be reasonably approximated 
by a classical one only in the large-S' limit BH . 




Figure 7: Critical temperature as a function of the decaying parameter p for a one- 
dimensional (a) and two-dimensional (b) spin- 1/2 Heisenberg ferromagnet. The dashed 
lines refer to the Green function (GF) method and to the modified spin-wave (SW) theory 
results, respectively. The dots represent MC predictions. 



From Eq. (I4.85t . we can easily obtain the behavior of the magnetization per spin 
m(r) = a(r)S, with T = T /JS^. For all values of p and d for which a phase transition 
occurs, we get 

fT -T\ 

m{T) ~ y—f^j , (4.95) 
implying (3 = 1/2 for the magnetization critical exponent. 
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The calculations are more delicate to obtain the near criticality behavior of the param- 
agnetic susceptibihty defined as x = 7% = J{<^/h) in the limit h with a ^ 0. In this 
case, the set of Eqs. (I4.58t -( l439l . with the appropriate expression for R(k) in Eq. (I4.60t . 
reduces to 



1 



1 



Pik) 



(4.96) 
(4.97) 



course, these equations can only be solved numerically. How- 
ever, an estimate of x{T) as T — > T^, can be obtained assuming, as usual, the low-k be- 
haviors ( RT^ for D.'^i'^k). Then, a solution of Eqs. K%i and ( ITWl . with x(T) ^ oo as 
r ^ r^, is found to be 



T-Tr 



(4.98) 



where the (/7,(3f)-dependent susceptibility critical exponent y{p,d) is given in Table I for 
the different regions of the {p — (i)-plane where LRO exists, together with other critical 
exponents to be obtained below. 



Table I 



Critical 
exponent 


3 J f 2t/ 

y<p<[d+2 


d<2 

d>2 




(/ < p < <j 


_ d + 2 


d<4- 
d>4 


2<d<4 


p>d + 2 
d = 4 


d>4 


l{p4) 


1/2 

p—d 
2d-p 




1/2 

lln 




1/2 
1 




1/2 

2 
d-2 


1/2 

lln 


1/2 
1 


5{p,d) 


P 

2d-p 
id-ip 
2d-p 




3ln 




3 




d+2 
d-2 


3ln 


3 


a{p,d) 




Oln 









d-4 
d-2 


Oln 






Table 1: Main (;7,^/)-dependent critical exponents in different domains of the (p — fif)-plane where LRO exists. The double vertical 
lines separate the regions where LRI's are active and SRI regime occurs. The symbol x/„ denotes the main (T — Tc)- or h- dependence 
with a logarithmic correction (for instance, x(r) ~ (T — Tc)~^ In [l/(T — Tc)], m{h) ^ h^/^ \\nh\^^^, etc.) 
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We now determine the behavior of the magnetization per spin m along the critical 
isotherm as /j — > 0. It is easy to see that, from the basic ME's (I4.58I) - (I4.60I) . the equation of 
the critical isotherm for small h and a assumes the form 

\-^o^-yT,Q,{o/h)=0 , (4.99) 

where ^ 

Qja/h) = (a/h) f ^rri , (4-100) 



and Pc{k) is obtained from Eq. (I4.97t with % replaced by a//i and T = Tc- Eqs. ( 14.991 ) 
and (14. loot can be solved numerically but a suitable estimate of a as /i ^ can be easily 
obtained working again in terms of the dominant contribution of (k) ask^O. Then, 
we get 

m = aSr^h^^^^P^'^\ (4.101) 

where the expressions of the critical exponent 5{p, d) in the different domains of the (p — d)- 
plane are given in Table I. 

Finally, we wish to calculate the critical exponent a{p,d) of the specific heat in zero 
magnetic field as T ^ for values of p and d such that a transition to the FM phase 
takes place. It is worth noting that, for this purpose, one cannot use the expression (14.701) 
for the reduced internal energy per spin. Indeed, it has a physical meaning only in the 
thermodynamic regimes where the decoupling {S^S^ ^) ^ (55) (5'^ can be used. A reliable 



calculation can be rather performed using the general expression (14. 14t with the decoupling 
procedure (I4.43t appropriate near the critical point, where o —^0, which takes into account 
the effects of the longitudinal spin correlations. From this expression, one can easily show 
that the reduced internal energy per spin near criticahty assumes the form 



'iBz {2nY (Or 



u^-ha--y(0)a^--Ta(\+a-) ,^ - , (4.102) 



where now is given by Eq. (14.581) with the appropriate R{k) in Eq. (I4.60t . With Eq. 
(14. 102t we have all the ingredients to calculate the zero-field specific heat near criticality. 
Indeed, working a.th—^0 and a — > with a/h = x{T) as T and taking into account 

Eq. ( 14.961) with P(k) evaluated as ^ — > 0, Eq. (I4.102t can be conveniently expressed, in 
terms of T, % and h, as 

lT-tx-\ym^t - ^(1 +y(o)x)x-i . (4.io3) 

Then, for the reduced zero-field specific heat Cj^^q{T) = Ci,=o{T)/S, we have 
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From this expression one can immediately determine the specific heat critical exponent 



OL{p,d). Indeed, with x(r) ~A(r 
I), Eq. (141041 yields 



Tc)-'^ orx(r) 



'A{T-T, 



In 



(r-r,) 



(see Table 



1 



1. ..^ = 



Oi„-l 



T-Tc) In 



1 



a{p,d) = l-y{p,d) <0, 



a = Oln . 



(4.105) 

From these relations and the exponent y{p,d) found before, one can immediately obtain 
the expUcit expression for a{p,d). The different values of a{p,d) varying p and d are also 
given in Table I. Other critical exponents can be obtained using the known scaling laws. 

The data collected in Table I show clearly some features which signal the degree of 
efficiency of the approximations made to lowest order in the CSDM. First, the critical ex- 
ponents coincide with those for the classical spherical model. Besides, within the region 
denoted with LRO (SRI) in Fig. 1, the critical exponents are exactly the same of those 
for the short-range Heisenberg model. This suggests that improving systematically the 
one-pole approximation for Aj((o) within the spirit of the SDM may produce better results 
beyond the mean-field approximation. 

As a conclusion we note that for J = 1 with 1 < p < 2, when the LRI's are active, the 
explicit expressions of the coefficients in the power laws (I4.95t . ( I4.98t . (I4.101t and (I4.105t . 
as functions of p, can be found in Ref. I.24J . 



B. Low-temperature paramagnetic susceptibility in absence of LRO For obtaining the 
low-temperature zero-field susceptibility in absence of LRO, namely for p>2d and d <2 

0. 



in Fig. 1, one can use again Eqs. (I4.96t -( l497t but bearing in mind that now 
Then, with the same procedure used before for determining the near critical behavior of 
the paramagnetic susceptibility as T ^ when LRO occurs, we easily find, for reduced 
susceptibility as T — > 0, 



~ < 



exp[A^A'73r] 

p-d 



d 
2-d 



p = 2d 

2d <p<2 + d 
p = 2 + d 
p>2 + d . 



(4.106) 



As usual, for d = \, all the coefficients can be explicitly calculated and Eq. (I4.106t becomes 



3 
47? 



exp [2%^/9T] 



(!) 



(p-i) 



3A(p) 



ifT^'ln 



(p-i) 



< 27 



Up-2)T- 



p = 2 

2<p<3 

p = 3 
p>3, 



(4.107) 
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where 



A{p) 



2(1)7 



9{p-\Y 



sin 



1 



-1 



1 



(4.108) 



The previous results for J = 1 appear consistent with the ones obtained in Refs. 1^ OTl 
for the corresponding quantum Heisenberg chain. It is worth noting an important aspect of 
the result in Eq. ( 14.1061 ) for p = 2d and hence in Eq. ( I4.107t for p = 2. This finding for 
reduced susceptibility, although consistent with the one obtained by means of the modified 
spin- wave theory {29} and the EMM for the two-time GF's L31J for quantum model, does 
not reproduce the well known exact Haldane result for = 1 l42ll which contains a factor 

1/2 

proportional to T . So, we are induced to argue that one must expect a <i-dependent 
power factor in T in Eq. (I4.106t which corrects the pure low temperature exponential 

1/2 

divergence reproducing T for J = 1. However, the lowest-order approximation in the 
CSDM, used in our calculations, has been able to capture the exponential divergence as 
r — > and, hopefully, one can improve systematically this result by means of higher order 
approximations. 



4.5 Moment Equations within the two 5-poles Approximation for the 
Transverse Spectral Density 



The next higher order approximation in the CSDM is to use the three exact ME's (I4.34t - 
(I4.36t and to choose a two 5-f unctions structure for the Aj((o) involving three unknown 
parameters (see Subsec. 3.2). A reliable possibility, suggested by quantum studies |4||5l0 

iHnnniEi, is 

Aj((o) = 27i|?i^^'5((0-(0j)+Xl^'5((0 + (0j)}, (4.109) 



so that, Eqs. (I4.34I) - (I4.35I) yield 



Mo{k)=2Nm, 



(4.110) 



where M^j(k) (v = 1,2) denote the moments in the right-hand side of the Eqs. (I4.35t and 



(I4.36t . respectively. These moments involve higher order correlation functions of the spin 
variables and hence, to close the set of ME's (I4.110t . it is necessary to use appropriate 
decoupling procedures. Below we suggest some of them which can be conveniently used 
for higher order calculations. 



4.5.1 Full Neglect of Longitudinal Spin Correlations 



Using in Eqs. (I4.35I) - (I4.36I) the Tyablikov-like decoupling procedure 



(4.111) 
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and the trivial ones (siSi) ^ (st) (si) and (siSiSi) ^ (^s) (5|) {si), in Eq. iZM 
we have 



M2(k) 



INmo^p''^ + ^ £ [/(P) - Jik - k'] 



2Nm ( (ol^-'^) 



+ 2^"^-l^[jik')-Jik-k') 



(4.112) 



(4.113) 



which close the set of ME's ( 14.1 101 ). In the previous equations d^^^ is given by Eq. ( 14.251 ) 
and use is made of the relation (14.371) which expresses exactly the transverse spin correlation 

function (^StS^Sj in terms of Aj((o). From Eqs. (I4.110t and (I4.112t - (l4.112t . we can 

conveniently write a set of two self-consistent integral equations for Mi (k) and M2 (k) and 
then determine the parameters coj and (v = 1,2) in terms of Mi(^) andM2(^). Indeed, 



by Eqs. ( 14. llOt . it is easy to show that 



and 



with 



and 



Mijk) 
2Nm 



{Ty) 



+j^j:\jik')-Jik-k') 



2Nm \ k J 



+ 2ci^P'^^'L[jik')-j(k-k' 



Mijk') 
M2(k')' 

Mi(k') 



M2{k') 



•4 



M2{k) 
2Nm ■ 



X^^) =Nm+{-\) 



iMi^f2Nm_\ 



\M2{k)) ' 



1,2). 



(4.114) 



(4.115) 



(4.116) 



(4.117) 



Of course, when Eqs. (14.1 14l) - (l4.115t are solved and hence Mi(k) and M2{k) are known, 
Eqs. ( I4.116l )-( r4. 1171 determine Aj((o), as given by Eq. ( 14. 1091 ). and all the thermodynamic 
properties follow from the general formulas established in Subsec. 4.1. In particular, for 
the dispersion relation coj of the undamped oscillations and the magnetization per spin, we 
find 



[Ty] 



and 



k k! 



mc^S > —. 

NiM2{k) 



M2(k') 



(4.118) 



(4.119) 
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We wish to underline that the previous results are physically meaningful only for {T ,h)- 
values in the nearly saturated regimes where, neglecting the longitudinal spin correlations, 
they constitute a reasonable approximation. 

Equations (14. 114t - (l4.119t show that it is convenient to work in terms of the single 
parameter 

Mi{k) 



x{k) 



satisfying the self-consistent equation 

x(k) 



t k V 



M2{k) 
j(k')-j(k-k')\ x{k') 



(4.120) 



J{k')-J{k-k') x{k') 



(4.121) 



A full solution of this equation is of course inaccessible and one must resort to numerical 
calculations. However, analytical results can be obtained, as usual, by expansion in powers 
of T with J ^g^d'^k' {...)/ {InY and appropriate constraints on the values of h 

k' 

assuring the convergency of the integrals (see discussion in Subsec. 4.4.1). Working, for 
instance, to first order in T, Eqs. (14.1 14t . (I4.115t and (14.1 19t yield 



{4.U2) 



which is exactly the same result (I4.63t obtained with the one 8-pole ansatz for Aj((o). Then, 
from Eq. (14.1211) . one easily has 



4k) - (oj \ 1 + 



%k^ 



9^l"\h)a^P>{k) + !F^''>{h,k) 



ip)i 



(4.123) 



where !Fy''\h,k) (v = 1,2) are defined by Eqs. (14.641) and (I4.65t and cOq^ = co^ 

Of course, an estimate of the integrals fFy^^ can be made using the ex- 
pansions (14.661) for Q.^P'>(k) as ^ — > along the lines outlined in Subsec. 4.4.1. 
Now, from Eqs. (I4.114t . (14. 115t and (I4.120t . a straightforward algebra yields 

M,{k)^2NSoy,^,l^l-^ ^^P\h) + ^(^f^[\h)Q.^P)Ck) + 0^i'\kk)) }, (4.124) 



^0,1 - 1^=0 



and 



M2{k)^2NS^ll^{\-i^_ 



T 



7i'\h) + ^ (^7i''\h)Q.^P)Ck) + fi''\h:k)) |. (4.125) 



Then, Eqs. (14.1161) and (14.1171) provide 

T 



%k^ 



)F[P\h)dP'^{k)+!F'f\h,k) 



(4.126) 
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and 

NS [1 + (-ir+i] 1 1 - L!F}P\h)'^ + 0{T^) 
lNm + 0{T^) , v = l 

0{T^) , v = 2. ^^-^^'^ 



The previous expressions for (Oj and X^^^ (v = 1,2) are the solution of the ME's (I4.114t in 
the nearly saturated regime to the first order in T by fully neglecting the longitudinal spin 
correlations. 

It is worth noting that, by inspection of Eqs. (I4.122t . (I4.126t and (I4.127t . the magne- 
tization and the dispersion relation, here obtained to the first order in T , coincide with the 
corresponding ones produced by the one-pole ansatz for Aj((o). Of course, possible differ- 
ences will occur beginning from second order in T. In any case, all the physical predictions 
of Subsec. 4.4.1 about the occurrence of LRO (see Fig. 1) remain qualitatively unchanged. 
All this is a clear signal of the internal consistency of the CSDM. 



4.5.2 Decoupling Procedures for Taking into Account the Longitudinal Spin 
Correlations in the Nearly Saturated Regime 

The previous approximation involving two 6-f unctions for Aj((o) can be systematically 
improved using in the ME's (14.1 101) a closure decoupling procedure which, although based 
on the near saturation relation S'j ~ S — SJ /2S, allows us to take properly into account 

the effect of the correlation function (^S^S^ k"}' perform this program, together with the 
Tyablikov-like decoupling (14.1 lit , for the higher order longitudinal correlation function in 
M2{k), we will use now the Hartree-Fock-like decoupling procedure 



CZ CZ CZ 



(4.128) 



With this prescription, the moment M2{k) in the right-hand side of Eq. ( I4.36t can be ap- 
proximated as 



M2ik) 



Ih^Nm - 4Nm^ 



JiO)-J{k) 



1 2 



+ 





4'\H'--k 


V 






k') 


|2(olp')+m 


7(0) 



(4.129) 
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At this stage, to close the system of ME's (I4.35t - (l4.36t . it is necessary to express the sums 
involving (^SiS'' in Mi(k) and M2{k) in terms of A^((o). Essentially, one must find a 
convenient approximate method to write sums of the type 



1 



(4.130) 



in terms of the original Aj((o). This is not a simple problem and here, within the spirit of 
the CSDM, we suggest a possible procedure which is suitable for exploring nearly saturated 
regimes. Assuming, ind eed, 5'^ . ~ S — SjSJ /2S and working in the Fourier space, with a 
simple algebra, the sum (14.1301 1 can be written as 



5j~iV%(0)5- 
1 



%(0) 



q 



+ 



452A^3 ^ 

qi,...,qi 



L %(-^3 - qA)^q,+q2;-q,-q, {S^S^.StS^^j , 



(4.131) 



which involves an higher order transverse correlation function. Now, with A = St S„^St^ 

^ yi q2 Hi 



and B = , we introduce the higher order transverse SD 

H,-rqM = i{{st(^)SqS^H(^)^Sq,} 

Then, in view of the general formalism developed in Sec. 2, we have exactly 



(4.132) 



+ 00 



^/(oA?,,...,54M 



271 



CO 



Hence, in Eq. (14.1311) . one can write, again exactly, 

12 ^qi+qi;-q?,-qA {^^^i^q^^X^q^ ^ = ^ 

where 



qi.qi 



cf(0 Ag3g^((0) 
271 CO 



A?3,?4(») = L hi+q2;~h~qAqu-AS^)- 



(4.133) 



(4.134) 



(4.135) 



qim 



As a next step, we perform in Eq. (14.1321 the natural and convenient decoupling procedure: 



A§.....a((0) ^ {slS-_r^)K-q. [K{^?3(^)'^?4} 



+ 



(4.136) 
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Then, the spectral function (I4.135t becomes 



^?3,$4(») 



?1 



13 



A-54(»)> 



(4.137) 



and, through Eq. (I4.134t . we have that the sum in Eq. (I4.131t can be written as 



1 



2^ L ^^(-^3-?4)8^,+§2;-§3-^4(5t5./t5.^^ ~ 



452A^3 



91 94 



1 



■I +%(^-/^)] (4^=9) H^--p 



(4.138) 



With this result and Eq. ( I4.131I ). it is easy to show that the general sum ( 14.1301 ). involving 
longitudinal correlation functions, can be expressed in terms of the transverse ones as 



Sl^Nm^O.^^{0) + {Sp--r) (sis 



9 -9/ \ p -p 



(4.139) 



where, in the Fourier representation, 



m - 5'- — ^ 



1 



5f=9 

25 



(4.140) 



Equation ( I4.139t is the main result for our next purposes. In terms of this expression and the 
exact relation (I4.37t for (s^S^j*"^ in terms of A-^(co), the moments Mi(k), with D.^{q) = 

J (q)- J (k-q), and M2ik), with aj^{q)= J{q)-j(k-q) ^2(£t^^ +m J{q)-j(k-q) |, 
can be easily expressed as 



^ ' k ' ^ '\.} -00 2% (St 

(4.141) 
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and 



M2{k) 



IKbTg) 



(Ty) 



N 



£ j(k')-j(k-k'] 



-oo 2% CO 



■ + 



+ 2^1 \jCki-k2)-jCk-Ch-k2)) 



ki.h 



{■■ 



J(k')-j(k-k') 



v=l 



271 CO 



(4.142) 



It is worth pointing that the previous ME's have been derived without assuming a particular 
functional structure for Aj(co). So, if we consider, for instance, the first two ME's (I4.34I) and 
( I4.35t with Ml (k) given by Eq. ( I4.14U and assume for A^(co) the one 5-pole ansatz (14.381 . 
it is possible to reformulate the analysis of Subsec. 4.4.1 but taking now into account the 
effects of longitudinal spin correlations. However, the aim of the last part of this subsection, 
as a suggestion for future more elaborate investigations, is only to show how the basic 
equations of previous subsection are modified using the two 6-poles ansatz (4.109) and the 
expressions (I4.141t and (14.1421) for M\ (k) and M2(k), respectively. 

First, we note that Eqs. (14.1 161) and (14. 117t . formally obtained by Eqs. (14. llOt . remain 
unchanged but now, for Mi(k) and M2{k), we must consider the system of self-consistent 
equations (to be compared with Eqs. (14.1 14t and (I4.115t ) 



My{k) 
INm 



and 



Mi{k') ^ 
M2(k') 



[J{k,-k2)-j{k-{k,-k2)) 



M,{k,)M,{k2) 

M2{h) M2(k2)'' 



(4.143) 



M2ik) 
INm 



+ [j{k,-k2)-J{k-{k,-k2)) 

k\k2 



Mjjk') ^ 
M2(k') 



2cop' +mJ{ki-k2)-Jik-(ki -h)) 



M,{h)My{k2) 

M2(ki)M2(k2) 



(4.144) 
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with m given again by Eq. ( I4.119t in terms of M\ (k) and M2(k). As we see, the appropriate 
variable for the self-consistent problem (I4.143l) - (l4.144t is again x(k) = Mi(k)/M2(k) for 
which it is immediate to obtain the appropriate equation. Then, when x(k) and hence Mi (k) 
and M2(k) will be known (analytically in some regimes, and numerically), the parameters 
which enter Aj((o) and m can be obtained using the Eqs. (I4.116t . (I4.117t and (14.1 lOt . In 
particular, for the dispersion relation we should have 



or 



(Ty) . 

K 



1 + i ^ £ - J(k - k')\x{k') + 



CO 



(Ty) N ' 



+ 



m 1 



52 / (7-v)\2A^2 



£ Jiki)-j(ki-k2) 



2(sLp^+m J(ki-k2)-j(k-(ki-k2)) x(^i)x(^2)} 



(4.145) 



k 



^+s^Ti]^2j;r^L[m-k2)-j{k-{ki-k2: 



k 



J{h-h)-j(k-{ki-k2))]]x(ki)x{k2)] 



(4.146) 



where we have denoted with co^'^^ the dispersion relation (14. 118t obtained with full de- 
coupling for the longitudinal spin correlation functions. All the remaining thermodynamic 
properties follow from the general formulas established before. 

Eqs. (I4.141t - (I4.146I) and, in particular, the Eqs. (I4.145t and (I4.146t for the frequency 
spectrum coj of the undamped elementary spin oscillations, show in a transparent way how 
the CSDM works taking into account systematically higher order effects on the relevant 
macroscopic quantities. This is just the key idea of the SDM. 



4.5.3 Suggestions for Future Studies 

The critical properties and the thermodynamic regimes with near-zero magnetization can 
also be explored within the two 5-functions representation (I4.109t for A^((o). This can 
be performed again on the basis of Eqs. (14.1 lOt with M2{k) given by the approximation 
(I4.129t which preserves the effect of the longitudinal correlations as the exact M\ (k) in Eq. 



(I4.35t . But now one must find an appropriate and convenient procedure to treat the sums 
in My(k) (v = 1,2) containing the correlation function ('^|'^^^/ for describing physical 
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situations where m > 0. A possibility consists in extending the decoupling procedure (I4.43I) 

as nnm 

«A^m2%(0)-i (l-^) ^L^ii^){sp-g)- (4.147) 

The next steps will be very similar to those made in Subsec. 4.4.2, but now much more com- 
plicated. This analysis is beyond the purpose of the present article. Our intention here was 
rather to give the main ingredients and suggestions for future possible analytical and numer- 
ical developments (also for nonmagnetic classical many-body systems), with the spirit of 
the SDM, for systematic higher order approximations satisfying more and more sum rules 
for the appropriate SD. In the next subsection, using some of the general results obtained 
above which are independent of the functional structure assumed for Aj((o), we will show 
how one can conveniently study the low-temperature damping of the spin oscillations along 
the lines suggested in Sec. 3. This will be performed using the general ME's (I4.34t - (l4.36t 
with Mv(^) (v = 1,2) in the form (I4.141t and (14. 142t for taking into account the effects of 
longitudinal spin correlations. 



4.6 Damping of the Spin Excitations in the Low-Temperature Regime Via 
the CSDM 

At this stage, we have all the ingredients to study the damping of the elementary spin 
oscillations by means of the extended CSDM suggested in Sec. 3. For our transverse 
SD Aj((o), according to the general prescription ( 13. 9t . we assume the modified Gaussian 
structure 



A^((o) = 27i(0?ir. 



exp 



(4.148) 



(7rr^)V2 

with the basic condition rj/tt>| ^ 1- This representation contains three unknown parame- 
ters ?i^,(Oj,r^ and hence we need to consider the three closed ME's 



L 2^^^(«)=2^-' 

y_ —(0A-,{(0)=M,{k), (4.149) 

^ —0i%{0i)=M2{k), 



where Mi (k) and M2{k) are assumed to be given by Eqs. (14. 14 It and (I4.142t . respectively, 
as obtained by means of the decoupling procedures (I4.111I) . (14.1281) and ( 14. 139t . They 
take into account the longitudinal spin correlation effects in the nearly saturated regimes in 
which we are interested here. 

Inserting the expression (14. 148t in Eqs. ( 14. 149t and performing standard integrals we 
have, for the parameters A,j,(Oj,rj, the equations 

X-^oy^ = 2Nm, (4.150) 
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+ 5^1 \J{h-k2)-j(k-{h-k2)) 



(4.151) 



4 + h 



2Nm((of-^y + J{k')-J{k-k' 



2 2(4^^) 



h + 



Here, m is given by 



+ I \j(k,-k2)-jCk-{h-k2)) 



2(ol^^ + m 



j(ki-k2)-J(k-(ki-k2)) 



k, k2 



(4.152) 



T ^ /^+" J(oAj((o) 

mC:iS— ^ 2_i 



ISN^h' J -oo 2k (0 



(4.153) 



and hence, with = 2Nm/(i>^ from Eq. (I4.150t . Eq. (I4.42t for the one 5-pole ansatz is 
reproduced. However, now co^ has a much more complicate structure. 

It is immediate to see that Eqs. (I4.151t and (I4.152t can be rewritten as (with 5W"v(^) = 
My(k)/{2Nm)) 



^k+2^kl^k = ^^^^) 



(0f + ^r^ = 3f2(^), 



(4.154) 



where 



Mi{k) 



k' 



j(k')-j{k-k') 



J{ki-k2)-J{k-{ki-k2)) 



ki ,k2 



(4.155) 
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and 



CO 



(Ty) 



J{k')-J{k-k') 



+ 



J{ky-k2)-J{k-{h-k2)) 



ki,k2 



k 



J{h-k2)-J{k-{h-k2)) 



(4.156) 



Now, consistently with the near saturation condition and with r^/co| <^ 1, for which the 
concept of damped oscillations has a physical meaning (see Subsec. 3.2), we show below 
that a solution of the self-consistent problem (14. 154t - (I4.156I) exists in the low-temperature 
regime. By simple manipulations of Eq. (14. 154t . we find 



M2{k)-M^{k), 



and, for /co? ^ 1 , we have formally 



T-^c^M2{k)-Ml{k). 
Thus, from Eqs. (14.1541 . it follows 

({&^'iM^(k)-lM2{k). 



(4.157) 



(4.158) 



(4.159) 



Of course, since F^ is expressed in terms of co^ through 9>{y,(k) (v = 1,2), our problem is 
reduced to solve the self-consistent Eq. (14. 159t for the dispersion relation co^ under physical 
conditions Fj > and T^/u^ <^ I. 

Focusing on the low-temperature limit and neglecting in Eqs. (14. 155t - (l4.156t terms 
> 3), with m\T=o = S and ro||T= 
(14. 1.5911 (with (0^ > 0) yield (as ^ oo) 

<t>(k;ki,k2) 



0{T",n > 3), with m\T=o = S and (o||r=o = [h + S{J{0) -j(k))f = ro^., Eqs. (I4l58l and 



IBZ 



{iTlY J (271 



IBZ 



h + S{J{0)-Jiki)) h + S{J{0)-Jik2)) 



(4.160) 



and 



^k- 



h + S{J{0)-J{k)) 7['\h)dP^(k)+7['\h,k) +0{T^), (4.161) 



where 



^(k;ki,k2) =J^\ \a'^P\k^-k2)-^'^''\k-(ki-k2)) 



+ 



D.^P\ki) - Q.^P\k -ki) a^P^(k2)-Q.^P^(k-k2) }. (4.162) 
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Of course, the solution (I4.160t - (l4.161t satisfies the conditions assumed at beginning. It is 
worth noting that, to the leading order in T, the frequency spectrum of the spin excitations 
coincides with the one obtained by means of the one 8-function ansatz for A^((o) (see Eq. 
(I4.62t ). A similar situation occurs for m for which the expression (I4.122t was found. Now, 
using the extension of the CSDM including the modified Gaussian functional structures 
for Aj((o), we have been able to determine the low-temperature damping of the classical 
oscillations. 

The integrals involved in the previous solution are very complicated for general d and 
p. An estimate could be obtained, as usual, using for j(k) or q!'P'> (k) the dominant contri- 
butions as 1^1 —>■ (see Eq. (I4.66t ). It is instructive, however, to present analytical results 
for case d = 1 and SRFs (p = oo) for which the integrals can be exactly performed. On finds 
indeed, with 7(0) = 27(1 -cosA:) {0<k< In) 



~ l6J^T^{\-coskf 



h{h + 4SJ) 



(1 -cos^) , 



""dk 1 
71 /j + 257(1 —cosk) 

2 



and 



(Ok ~/i + 475'sin^ ( - 



1 



4752 



1 



h 



h + ASJ 



(4.163) 



(4.164) 



Notice that, due to the approximation used for S] under nearly saturation conditions, the 
previous results break down at any given low temperature where the external magnetic field 
becomes sufficiently small (see discussion in Subsec. 4.4.1). In particular, for J = 1 and 
/? = 00, they have a physical meaning only ilT ^h^/^ for very low values of h (see Eq. 
dH). 

It is interesting to compare the results (14. 163t and ( I4.164I) . here obtained by means 
of the extended CSDM, with those obtained by different methods applied directly to the 

Heisenberg ferromagnetic chain with short-range interactions. First, the Gaussian peak 

1/2 

width exhibits the correct low-temperature behavior (oc T) found by other authors 
l40l l43l . Furthermore, our result for (Ok is consistent with the corresponding exact one 
obtained in Ref. |43l. Of course, further comparisons for the cases d > \ and p ^ °° 

could be performed by means of numerical calculation taking also into account that we find 

1/2 

r4 oc T for any d and p with h ^ 0. 



5. Conclusions 



In this article we began with a review of the general formalism of the two-time GF's tech- 
niques in classical statistical physics in a form quite parallel to the well known version 
in quantum many-body theory. Particular emphasis has been given to the spectral prop- 
erties and to the concept of SD which plays a relevant role for our purposes. The EMM 
for calculation of classical two-time GF's is also presented but we have mainly focused on 
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the general ingredients of the less known classical version of the SDM. This method, intro- 
duced by Kalashnikov and Fradkin L4 5 1 many years ago in quantum statistical physics, has 
been extensively tested by several authors for a wide variety of quantum many-body sys- 
tems |4 3'.0|8l|9|[ro|[n- We have explained in detail, also with a pedagogical aim, how 
the CSDM works in classical many-body physics for obtaining systematic nonperturbative 
approximations to explore the macroscopic properties of a wide variety of systems, also 
including phase transitions and critical phenomena, without the explicit calculation of the 
(generalized) partition function. Besides, general insights and key ideas, for an implemen- 
tation of the method to investigate the dispersion relation and the damping of the classical 
excitations on the same footing, are also given. The possibility to formulate a many-body 
theory in the framework of the classical statistical mechanics offers the remarkable advan- 
tage to draw insights directly from the wide baggage of powerful mathematical techniques 
and approximation methods developed in many years of research activity in quantum many- 
body theory which have sensibly contributed to the great progress in condensed matter 
physics and, more generally, in quantum statistical physics. It is also worth noting that the 
classical version of both the EMM and the SDM, with the use of commuting dynamical 
variables, allows to avoid the intrinsic difficulties related to the noncommutability of the 
operators which are the basic elements of the quantum framework. This constitutes a very 
favorable opportunity to have physically relevant information about the macroscopic prop- 
erties of microscopic models when the quantum fluctuations do not play an important role. 
Of course, the use of classical dynamical variables and the machinery of a reliable classical 
many-body framework is also particularly advantageous for numerical simulations. 

Next, the effectiveness of the CSDM has been explicitly checked with the study of the 
thermodynamic properties of a highly nontrivial classical J-dimensional spin-S' Heisenberg 
FM model with LRFs decaying as r^P (p > d) in the presence of a magnetic field. However, 
the method is general and may be conveniently used for studying other classical many -body 
systems especially when a reliable nonperturbative method is required. Focusing on the 
selected spin model and working in some detail to lowest order in the CSDM, a rich physical 
scenario has been obtained. The main magnetic quantities in the nearly saturated regime 
have been calculated analytically as functions of d and p. The FM LRO at finite temperature 
has been shown, in a transparent way, to take place in a wide region of the (p — <i) -plane 
consistently with exact and numerical predictions. The thermodynamic regimes with zero 
and nearly zero magnetization have been also investigated by means of a reliable decoupling 
procedure, as suggested by quantum studies |5 37|, which takes into account the effect of 
longitudinal spin correlations. This allowed us to estimate the critical temperature, the 
critical properties in terms of d and p beyond the mean field approximation and the low- 
temperature behavior of the paramagnetic susceptibility when LRO is absent. 

All this supports our conviction that the CSDM, as in the quantum case, constitutes 
a very effective method for describing the macroscopic properties of classical many-body 
systems. It is indeed surprising that, already to the lowest order of approximation, the 
method is able to capture, in a relatively simple way, the relevant physics of an highly 
untrivial spin model, consistently with MC simulations 1211 l22l l32l and exact results IT6l 
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[T7l[T8l[T9l|20l|36l|42l|. It is also worth pointing that the method provides the possibility 
to perform, systematically and in a transparent way, higher-order approximations with the 
aim to improve the results here presented. Some suggestions along this direction, also for 
exploring the damping of the classical excitations, have been given in Sees. 4.5 and 4.6, but 
new ideas can be profitably drawn from the big experience acquired working by means of 
the quantum version of the SDM. 

In conclusion, we believe that the CSDM and the formalism of the two-time GF's in 
classical statistical mechanics constitute a powerful tool to study efficaciously the equilib- 
rium and transport properties of classical many-body systems, also in nonextensivity condi- 
tions [15 1 . Further developments and applications are, of course, desirable and we warmly 
hope that the present article will stimulate further and more elaborate and meaningful stud- 
ies. 
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